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Ideal triangle groups, dented tori, 
and numerical analysis 

By Richard Evan Schwartz* 



Abstract 

We prove the Goldman-Parker Conjecture: A complex hyperbolic ideal 
triangle group is discretely embedded in PU(2, 1) if and only if the product of 
its three standard generators is not elliptic. We also prove that such a group 
is indiscrete if the product of its three standard generators is elliptic. A novel 
feature of this paper is that it uses a rigorous computer assisted proof to deal 
with difficult geometric estimates. 



1. Introduction 

A basic problem in geometry and representation theory is the deformation 
problem. Suppose that po : T ^ Gi is a discrete embedding of a finitely 
generated group T into a Lie group Gi. Suppose also that Gi C G2, where 
G2 is a larger Lie group. The deformation problem amounts to finding and 
studying discrete embeddings : F — > G2, which extend po- 

When r is the fundamental group of a surface, Gi = Is(H^), the isom- 
etry group of the hyperbolic plane, and G2 = Is(H^), the isometry group of 
hyperbolic three space, one is dealing with the theory of quasifuchsian groups, 
which is quite well developed. (See, for instance, the bibliography in [T].) 

The complex hyperbolic plane CH^ is a 2-complex-dimensional manifold, 
which is negatively curved, Kahler, and also a symmetric space. It contains 
as a totally real, totally geodesic subspace, and is often considered to be its 
complexification. The theory of deforming Is(H^) representations into Is(CH^) 
is also quite rich, though much less developed. (For a representative sample of 
such work, see [FZ], [GKL], [GuP], [KR], [Tol].) 
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1.1. The Goldman- Parker conjecture. In [GP], Goldman and Parker took 
one of the first steps on the road to a theory of complex hyperbolic quasifuch- 
sian groups. They defined and partially classified the complex hyperbolic ideal 
triangle groups, which are representations ps : F — > Is(CH^). Here F is the free 
product Z/2*Z/2*Z/2. The representation ps maps the standard generators 
of F to distinct, order-two, complex reflections, such that any product of two 
distinct generators is parabolic. A complex reflection is an isometry which has, 
for a fixed point set, a complex line in CH^. (See §2.3 for more information.) 

Modulo conjugation, there is a one-parameter family {ps\ s G R} of such 
representations. The indexing parameter, s, is the tangent of the angular 
invariant of the ideal triangle formed by the three complex lines fixed by the 
generators. The angular invariant measures the extent to which the vertices 
of the triangle fail to lie in a totally real subspace. (See §2.5 for a definition.) 
The representation po is the complexification the familiar real ideal triangle 
group generated by reflections in the sides of an ideal triangle in H^. 

Let s = (105/3)^2 and s = (125/3)^^. According to Goldman-Parker: 

Theorem [GP]. // |s| > s then ps is not a discrete embedding. If\s\ < s 
then Ps is a discrete embedding. 

Here is a sketch of the first half of this result. Let gg be the product of 
all three generators of Ps(X), taken in any order. In [GP] it is shown that 
is loxodromic for |s| G [0, s), parabolic for |s| =s, and elliptic for |s| > s. (See 
§2.1 for a classification of isometrics.) If gg is elliptic, with finite order, then ps 
is not an embedding. If gg is elliptic, with infinite order, then ps is not discrete. 

The theorem above is not sharp. For jsj > s the analysis in [GP] breaks 
down, but Goldman and Parker conjecture that ps remains a discrete embed- 
ding for \s\ G {s,s]. The significance of the result /conjecture combination is 
that it proposes the first complete description of a complex hyperbolic defor- 
mation problem. 

1.2. Results and methods. The purpose of this paper is to prove a sharp 
version of the Goldman-Parker conjecture. 

Main Theorem. ps is a discrete embedding if and only if gg is not 
elliptic. Also, ps is indiscrete if gg is elliptic. 

Our indiscreteness proof uses some Galois theory to show that gg must 
have infinite order when it is elliptic. 

For our discreteness proof, let Tg = ps(^)- Let 5CH^ be the ideal bound- 
ary of CH^. We construct a surface-like set Z(s) C dCH^, which we call a 
dented torus. We then prove that the orbit TgZ{s) consists of (essentially) 
disjoint surface-like sets. This phenomenon feeds into a variant of the Klein 
combination theorem to prove that pg is a discrete embedding. 
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We found Z(s) experimentally, using our interactive program [SI]. The 
technical part of this paper consists in making the pictures in [SI] into rigorous 
proofs, via computation and numerical analysis. While our proof is logically 
independent from [SI], most of the ideas came from analyzing the output of 
this program. 

We use two kinds of programs to verify our computational claims. The 
first program, which computes quantities depending on the eigenvalues of gg, 
is Mathematica code [M]. This code runs with much more precision than we 
need. The second program, which is used for the bulk of the computation, 
is C code [KeR]. Our C code uses only those operations which conform to 
the IEEE standards [I] and employs interval arithmetic to bound the roundoff 
errors rigorously. Our code takes about 350 computational hours, when run 
on a Sparc Ultra 5. We ran the code, in parallel, on 11 such machines, over 
the course of a weekend. 

Our code includes a variety of interactive visual displays, which allow the 
user to see, to a large extent, that everything operates as claimed. The code 
is available upon request. 

We give explicit parametrizations of the main objects in the paper, as well 
as pseudocode expositions of the main algorithms. We hope this is enough 
information for the inclined (and energetic) reader to reproduce the results 
here in an independent experiment. 

1.3. Overview. In Section 2 we review some basic complex hyperbolic 
geometry. 

In Section 3 we define the triangle groups, establish the indiscreteness part 
of our main theorem, and study the eigenvalues of gs- 

In Section 4, we study the action of Fg on the Clifford torus, 

T = {{z,w) e C^l \z\ = \w\ = V2/2} C S^. 

(The ideal boundary of CH^ can be naturally identified with the three-sphere.) 
We give an alternate proof of the Goldman-Parker discreteness result. Though 
our proof breaks down at the same place the proof in [GP] breaks down, it 
sheds new light on what should be done in the critical interval (s, s]. Our basic 
idea is to put some dents in the Clifford torus, so that the resulting surface 
moves around better under the action of the group. 

In Sections 5 and 6 we construct the "surfaces" which replace the Clifford 
torus. In Section 5, we introduce a projectively natural coning operation, the 
hybrid cone construction. This procedure creates surfaces called hybrid sectors, 
which are certain embedded solid triangles. In Section 6 we introduce, for each 
parameter s, the dented torus, Z(s), a surface which is made from cutting out 
two solid triangles from the Clifford torus and gluing back six hybrid sectors. 
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In Sections 7 and 9 we give the bulk of the proofs. In Section 7 we give 
the discreteness proof of the Main Theorem for the parameter s, modulo two 
technical estimates, 1 and 2. In Section 8 we establish the Year Lemma, a 
result which shows that we only need consider 365 parameters in the interval 
J = [s, s] . In Section 9 we give the general discreteness proof, modulo the two 
estimates, 1 and 2. 

The remainder of the paper is devoted to establishing and illustrating 
Estimates 1 and 2. In Section 10 we describe the method by which we estimate 
the location of a hybrid sector in space using finitely many computations. In 
Section 11 we detail the computation scheme used to verify Estimates 1 and 
2. In Section 12 we discuss the implementation and computational accuracy 
of our code. Wc also give a record of the actual calculations. In Section 13 we 
illustrate Estimates 1 and 2 with computer plots, produced with our code. 

Suggested reading. We have indicated, through a series of remarks, the 
material in Sections 2-6 which is strictly necessary for the proof of the parabolic 
case. Read only this material, and then read Section 7. Following this, skip to 
Section 13, and look at the computer pictures, which give compelling evidence 
that Estimates 1 and 2 are true. Afterwards, try to understand the general 
case. 

1.4. Related results. This paper is the first in a series. We prove in [S2] 
that the orbifold at infinity, in the parabolic case, is commensurable to the 
Whitehead link complement. [S2] also gives a noncomputational proof that 
there is some e > such that ps is a discrete embedding for s G (s — e, s]. 

In [S3] we explain (among other things) how to build models out of string, 
of the limit sets of the complex hyperbolic ideal triangle groups. 

In [S4] wc consider an extremely deformed discrete representation of the 
(4, 4, 4)-Hecke triangle group, and use it to produce the first example of a closed 
hyperbolic 3-manifold which bounds a complex hyperbolic 4-manifold. 

For some neat results and conjectures on the algebraic side of the complex 
hyperbolic ideal triangle groups see [Sa]. 

For deformations of other triangle groups into PU(2, 1), see [W-G]. 

For another example of an application of computation to hyperbolic ge- 
ometry, which is independent from ours, see the monumental [GMT]. 

1.5. Acknowledgements. I must say that writing this paper has been an 
extremely gruelling experience, though sometimes an exhilarating one. I would 
like to thank Bill Goldman for his perspective and encouragement as well as 
for the efficient proof of Lemma 2.2. I would like to thank Martin Bridgeman, 
Peter Doyle, Jeremy Kahn, Josh Maher, Andy Mayer, Robert Miner, Robert 
Meyerhoff, John Parker, Bill Thurston, and Justin Wyss-Gallifent for inter- 
esting conversations relating to this work. I would like to thank the tireless 
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referees for many helpful stylistic suggestions. Finally, I would like to thank 
my wife Brienne Brown for her C programming advice, and in particular for 
suggesting the elegant subdivision algorithms used in Section 10.8 and Sec- 
tion 11.1. 

I dedicate this paper to my daughter, Lucina Caroline Schwartz, who 
was born right around the time I started thinking about triangle groups, and 
who is now learning about triangles herself. I also dedicate this paper to 
the memory of Hanna Sandler, a great enthusiast for the complex hyperbolic 
triangle groups, who passed away this winter. 

2. Background 

2.1. The complex hyperbolic plane. C^'-*^ is a copy of the vector space 
equipped with the Hermitian form 

(1) {U, V) = UiVi + U2V2 - U3VS. 

The spaces CH^ and 5CH^ are respectively the projective images, in the 
complex projective plane CP^, of 

(2) N- = {veC^'^\ {v,v) <(}}; No = {veC^''^\ {v,v) = 0} 
(See [G, p.67] or [E].) The map 

(3) @ivi,V2,V3) = (—, — ) 

\V3 V3J 

takes and Nq respectively to the open unit ball and unit sphere in C^. 

Henceforth we identify CH^ with the open unit ball. 

Given a point V G CH^ U S'^, we will say that V £ Q~^{V) is a lift of 
V and that a lift of the form (fi,f2, 1) is affinely normalized. We define the 
vector Q~^(y) as the afHnely normalized lift of V. 

SU(2, 1) is the group of (, ) preserving, determinant 1 complex linear 
transformations. PU(2, 1) is the projcctivization of SU(2, 1), and elements of 
PU(2, 1) preserve CH^. Concretely, each T G SU(2, 1) determines an element 
of PU(2, 1) via 

(4) r = eofoe~\ 

The map T^TisaS — 1 surjective lie group homomorphism. 

An element T G PU(2, 1) is called loxodromic if T has exactly two fixed 
points in S'^ , parabolic if it has exactly one fixed point in 5^, and elliptic if it 
has a fixed point in CH^. This classification is exhaustive and exclusive. T 
is called ellipto-parabolic if T is parabolic and also stabilizes a complex line 
in CP^. (See [G, p. 203] for more details.) For instance, the element gj is 
ellipto-parabolic. 
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CH^ has a PU(2, l)-invariant Riemannian metric, unique up to scale. Let 
S denote the path metric induced by this Riemannian metric. When normahzed 
as in [G, p. 78], 

,2^iv,w) {V,W){W,V) - 

(5) cosh — - — = oiv, w). 
^ ' 2 {v,v){w,w) ^ ' ' 

2.2. Heisenberg space. Let 3^(z) and '^(z) denote the real and imaginary 
parts of a complex number z. The Siegel domain is the affine subset 

(6) {{z, w)\ SRH > C C CP2. 

Let H stand for C x R. We call H Heisenberg space. Given p e S^, a 

Heisenberg stereographic projection from p is a map B : 5^ — {p} H which 
has the form B = tto/?. Here tt{z, w) = (z, and /3 is a complex projective 

transformation of CP^ which identifies CH^ with the Siegel domain. We write 
GO = B(p) in this case. For instance 

m B„(.,») = (^,9l^) 

is a Heisenberg stereographic projection which maps (—1,0) to oo. The map 
Bo is a close relative of the Cayley Transform, given in [G, p. 112]. The map 
Bo conjugates the PU(2, 1) stabilizer of (—1,0) to certain affine maps of Ti. 

We call £"0 = {0} X R the center of Heisenberg space. 

The complex lines tangent to form a canonical contact distribution on 
S^. We call this distribution £. Heisenberg stereographic projection maps £ 
to a corresponding distribution in TC, which we give the same name. In Ti., the 
distribution £ is the null distribution to the contact form 

(8) (J = 2ydx - 2xdy + dt 

(see [G, p. 124]). In this formula, C has been identified with R^ in the usual 
way. £ has cylindrical symmetry. It is preserved under maps of the form 

{z, t) {\z, \\\h + r); A G C*, r G R. 

If n is a plane in £, based at the point {z,t), then the maximum slope of a 
vector in li is 2\z\. 

We say that a curve, in either 5*^ or H., is £-integral if its tangent vector 
at every point is contained in £. The term £-integral is sometimes called CR 
horizontal. We prefer to use the former term, because the term horizontal 
comes up elsewhere in the paper. 

2.3. C-circles. A complex slice is the intersection of a complex line in 
with CH^. Complex slices are totally geodesic subspaces, when considered 

as Riemannian subspaces of CH^. A C- circle (also known as a chain) is the 
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intersection of a complex line with S^, provided this intersection is more than 
a single point. A C-circle is a round circle, transverse to 6. A C-arc is a 
nontrivial arc of a C-circle. Given two points p ^ q E S^, there is a unique 
C-circle containing p and q. 

Here is a method for parametrizing C-arcs. Suppose X, Y, Z are three 
distinct points contained in a common C-circle, C. After switching the order 
of X and Y if necessary, there are lifts X,Y,X, unique up to multiplying all 
three vectors by the same unit complex number, such that 

(9) {X,Y) = i; {X,Z) = {Y,Z) = 1 
(compare §2.5 below). Note that iX + tZ e Nq for t G R, so that 

(10) Cz{X,Y;t) = eiiX + tZ), te[0,l], 

parametrizes the C-arc of C which joins X to Y and which avoids Z. 

Let = C^'-*^ — iV_ — Nq. If C is a C-circle, then there is the polar vector 
C* G Ar+, unique up to scaling, such that C = {v e No\ {v,C*) = 0}. There is 
a unique involution Ic G PU(2, 1) fixing C. This map is computed by setting 
Ic = @ o Ic* ° 0~^, where 

(See [G, p. 70].) Another feature of the polar vector is this: Two C-circles, Ci 
and C2 are linked in if and only if S{Ci, C|) < 1, where 6 is as in Equation 5. 

We say that a Heisenberg chain is the image of a chain under a Heisenberg 
stereographic projection. Curves of the form {{z} x R) U 00 are chains. In 
particular, the center of 7-^ is a chain (with 00 deleted.) All other Heisenberg 
chains are ellipses which project to round circles under the projection ttc : 

C. (See [G, p. 125].) Let C be such a chain, with center of mass m. Let 
Em be the plane of £ based at m. Let E„i be the affine plane which is tangent 
to Em at m. We call Em the prolongation of Em- We have C C Em- We will 
say that a round Heisenberg chain is one which is, itself, a round circle. The 
center of mass of a round Heisenberg chain is contained in the center of H- 
Such curves are contained in planes of the form C x {t}. 

2.4. Ti-circles. A real slice is a totally real, totally geodesic subspace of 
CH^. Every real slice is isometric to the real slice R^nCH^. An R-circle is the 
accumulation set, in S^, of a real slice. Every R-circle is PU(2, l)-equivalent 
to the particular R-circlc R^ n S^- All R-circles are £^-integral. 

We say that an R-circlc is round if it is a round circle. An R-circlc is 
round if and only if the real slice containing it lies in a real 2-plane in C^, 
which is true if and only if the real slice contains (0,0). The geodesies in a 
round R-circle are line segments in C^. Not all R-circles are round. An H-arc 
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is a nontrivial arc of an R-circle. There is an S* -family of R-arcs containing 
two points p q e S^. The union of ah such arcs is cahed a spinal sphere. 
(See [G, p.152].) 

Here is a method for parametrizing R-arcs. Suppose X, Y, Z are three 
distinct points contained in a common R-circle, R. There are lifts X,Y,X, 
unique up to scaling all three vectors by the same complex number, such that 

(12) {X,Y) = {Y,Z) = {Z,X)eIL. 

(Compare §2.5 below.) We call {X, Y, Z} an equalized choice of lifts. The 
curve 

(13) Rz{X,Y;u) = G{{l-u)X + uY + u{u-l)Z), uG[0,1], 

parametrizes the R-arc of R which joins X to Y, but which avoids Z. 

We say that a Heisenberg R-circle is the image of an R-circle under a 
Heisenberg stereographic projection. Any Heisenberg R-circle 7 which contains 
00 is (the extension of) a straight line. We call these R-circles straight. Note 
that 7 has the form (L x {t}) U 00 where L G C is a line through the origin, if 
and only if 7 is straight and intersects the center of 7^. In this case we call 7 
level. All other Heisenberg R-circlcs arc curves, which project to lemniscates 
via the projection vrc- (See [G, p. 139].) 

Here is a geometric interpretation of our parametrization: 

Lemma 2.1. Let Bq be as in Equation 7. If Z = (—1,0) and X,Y,Z lie 
on a common R-circle then Bo{Rz{X,Y;u)) = (1 — it)Bo(X) -|- uBo{Y). 

Proof. Since Bq maps the R-circle containing X, y, X to a straight line, 
the point P = 'Bo{Rz{X,Y;u)) must be a convex combination of Bo(X) and 
Rq(Y). Our parametrization of R-arcs is natural in the following sense: If, for 
j = 1,2, the points Xj,Yj,Zj are contained on an R-circle, then an element 
of PU(2, 1) which carries {Xi,Yi, Zi}, in order, to {X2,Y2, Z2}, also carries 
Rzi{Xi,Yi;u) to Rz.2{X2,y2;u). Also, Bq conjugates the stabilizer of Z = 
(—1,0) to affine maps of The facts just mentioned imply that the particular 
convex combination depends only on u, and not on the triple of points chosen. 

We will make the computation for the points X = (1,0), Y = (0, 1) and 
Z = (-1,0); an equalized choice of lifts is X = (1,0,1), Y = (0,2,2) and 
Z = (—1,0, 1). Equation 13 gives 

Rz{X,Y-u) = (\^,-^). 



We have Bo(X) = (0,0) and ^o{Y) = (1,0) and, from the calculation above, 
Bo(P) = u. Since the result is true in this case, it is true in all cases. □ 



IDEAL TRIANGLE GROUPS 



541 



2.5. The angular invariant. See [G, pp. 210-214] for proofs of the results 
in this section. 

Given three points X ^ Z E S^, the Cartan angular invariant is 

(14) A{X,Y,Z) = ^Tg[-{X,Y){Y,Z){Z,X)] G [-|,|]. 

We have A(X,Y, Z) = ±7r/2 if and only if these points lie in a common C- 
circle, and A(X, Y,Z) = if and only if these points lie in a common R-circle. 

2.6. The group law on S^. Recall that is a Lie group. With as the 
unit sphere in C^, as above, the group law on is 

(15) {zi,Wl) ■ iz2,W2) = {ziWi - Z2W2, Z1W2 + Z2W1). 

The right multiplication map RxiY) = Y ■ X \s an element of PU(2, 1) which 
acts isometrically on S^, in the round metric, and moves all points by the same 
amount. The point here is that Rx commutes with the complex structure, 
which is left multiplication by i. 

2.7. The Clifford torus. We do not know a reference for the material in 
this section. 

The Clifford torus is the torus T = {{z,w)\ \z\ = \w\ = V2/2} C S^. 
Define arg : T — Tq by the formula arg(2;, u;) = (argz, argu;). The preimages, 
under arg, of curves of slope 0, 1 and cxd are all chains on T. We call such 
chains, respectively, horizontal, vertical, and diagonal. The preimages, 
under arg, of curves of slope —1 are round R-circles on T, having the form 
{{z,w)\ zw = X}nS^, where |A| = 1. 

We say that a Clifford triangle is a solid triangle in T whose edges are 
contained in chains, one from each of the three foliations. See Figure 3.1 for 
an example. 

We say that a Clifford rectangle is a solid rectangle whose boundary con- 
sists of two vertical C-arcs and two horizontal C-arcs. In parametrizing a 
Clifford rectangle we wish to avoid the exponential function, for computa- 
tional reasons. Here is how we do it. Let Q be a Clifford rectangle. Let H and 
V be maximal horizontal and vertical C-arcs in dQ. 

Let Hi and H2 be the endpoints of H and let Hs be the midpoint of 
H — H. Here H is the C-circle containing H. Likewise define Vi, V2 and V3. 
We parametrize Q by 

(16) Q{t,u) = i7ri{CHsiHi,H2;t)),7r2{CvsiVi,V2;u))). 

Here tti and 7r2 are projections onto the first and second factors, and the 
components of our parametrization come from Equation 10. 
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V3 

Figure 2.7 

We call this the rational parametrization of Q. It appears that there 
are two choices for H and two for V. However, recall from our definition of 
CnsiHi, H2;t) that the points Hi, H2, must be correctly ordered. Like- 
wise for Vi, V2, V3. These extra constraints make the rational parametrization 
unique. 

2.8. Pure translation flows. 

Remark. The reader interested only in the parabolic case of the main 
theorem may skip this section. 

Say that a pure translation along a geodesic 7 G CH^ is one which pre- 
serves all real slices containing 7. Say that a pure translation flow along 7 is 
a flow (pt '■ CH^ CH^ such that cpt is a pure translation of 7, for all t. Let 
/ be the geodesic (—1, 1) x {0}. Let || || be the Euclidean norm on C^. 

Lemma 2.2. The vector field V{z,w) = (1 — z'^,zw) generates a pure 
translation flow along I. 



Proof. Let Bq be the map defined in Equation 7. Let (C)^^) = ^o{z,w). 
We compute 

v2 



dBo{V) 



-2 



1 

1+z 





z- -1 
zw 



w 



2 - 2z, 



z + V 1 + z 



(C,2a;). 



This vector field generates the flow iIj^{C,,uj) = {rC,,r'^uj). Clearly, is a flow 
by complex linear maps which preserves the Sicgcl domain. Hence, (pt is a flow 
by complex hyperbolic isometrics, obviously preserving / and R^. □ 



Lemma 2.3. Suppose that R is a real 2-plane which contains (0,0). Sup- 
pose pi,P2 £ R — CH^ — S*^ are two points such that (0, 0),pi,p2 o'^e all con- 
tained on a line L C R. Let (3 he the pure translation along L n CH^ such that 
P{Pi) = P2- Then, for all q 1^ S^, 



\\Pl -P2\ 



S = inf(||pi||, \\p2\ 
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Proof. The stabilizer subgroup of (0, 0) in PU(2, 1) acts isometrically on 
and transitively permutes the real 2-planes through (0,0), as well as the 
real lines through (0, 0). Thus, we may assume that ^ = R^, that L = I, that 

Pi = ixi,0) and p2 = (a:2, 0), with 1 < xi < X2 < oo. Let V be the vector field 
from Lemma 2.2. Since sup^-g^s ||y(x)|| = 1, we conclude that p{q,4't{Q)) < t 
for all q £ S^. The value of t, for which = P2, is 

_ /"^a dx X2 - xi _ IIpi -P2II 

X2 - 1 - xf - 1 " 51 - 1 • 



3. The triangle groups 

3.1. Basic definitions. Given s G [0, 00) let be the Clifford triangle 
having vertices 

(s + i) 



(17) {(3s, Ps); {13s, Ps); {f3s,(3sy, Ps 



V2 + 2s2 ■ 



The angular invariant of these three points is iarctans, depending on the 
order in which they are taken. The first part of Figure 3.1 shows the images 
of these points under the arg map. 




(0,0) 

Figure 3.1 



Let -ff+,s, V-\-,s and D be the horizontal, vertical, and diagonal C-circles 
containing pairs of these points. Let R be the R-circle which bisects A^. Define 

(18) d{z,w) = {w., z); r{z,w) = {w,z). 

Note that d is the complex reflection fixing D and r is an antiholomorphic 
involution fixing R. Note that d and r commute. 

We define H^ g = d{r{H^^s)) and V-^g = d{r{V+^s))- The second part of 
Figure 3.1 shows these chains. Let h±^s be the complex reflection fixing H± s, 
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etc. The group Tg = ps(^) is generated by the triple {d, hj^^g-, These are 

not the standard generators of F^, but they are quite useful to us. We define 

(19) Qs = /i+,s o h_^s ° d. 
Since r conjugates h±^s to v±^s and vice versa, 

(20) rT,r-^ = Ts; rg,r-^ = g-\ 

The representations ps and p-s have the same image in PU(2,1), so we are 
justified in taking s G [0, oo). We define J = [s, s]. Usually we take s & J. 

3.2. Matrices. We now give formulas for SU(2, 1) lifts of some of the 
elements defined above. 



(21) 









3 

4/3. 





-3 









3 





-4/3. 
-3 





■ 


-1 


■ 


(22) d = 


-1 










. 





-1. 


Setting g = h-^.h-d, we have 









(23) 





-Ms) 

-A2{s) 










Ms) 

-A^{s)\ 



M{s) 



s + 17i _ 
s + 1 



M{s) 



12^/2i 



Clearly, d, is a lift of the map d{z, w) = (w, z). One can check easily that 
h±^s is a determinant-one involution which preserves Nq C C^'"*^. One can also 
check that this matrix has, as an eigenvector, a vector which is polar to H±^s- 
One can also take traces of these matrices and their products, checking that 
they agree with the traces of the corresponding matrices in [GP] . 



3.3. Indiscreteness proof. We now show that p. is indiscrete for s G (s, oo). 
The work in [GP] shows that gg is elliptic for such s. We will show that gg has 
infinite order. 

Prom the formula for gg , we see that the trace of gg satisfies the equation 

(24) ^rg + N + l){Tg + N + l) = N^; Tg ^ -1 

for N = 8. Note that 3?(ts) < — 1 as long as G N. If gg has finite order then 

(25) Tg = uj^ + LJ^t + uj'n, = cxp(27ri/n); p + q + r = 0. 

Here n is taken as small as possible. Obviously n > 2. We suppose that 
Equations 24 and 25 are true for some N eN, and derive a contradiction. We 
set r = Tg. 
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Lemma 3.1. n divides 12. 

Proof. For k relatively prime to n, let (T^ : Q[ti;n] Ql^^n] be the Galois 
automorphism determined by cjfc(ti;„) = w^. Equation 24 is defined in Q[a;„], 
so it remains true if we replace r by crfc(r). In particular, 3fi((Tjk(r)) < — 1. 
Summing over all those k G {1, . . . , n — 1} which are relatively prime to n, we 
get 3fi(X]fc <7fc(T)) < —<p{n). Here (/> is the Euler phi function. Hence 

(26) |^c7fe(r)|></.(n). 

k 

For m = p,q,r let (m, n) be the greatest common divisor of m and n. Let 
c^m = n/{m,n). Note that wj^ is a primitive c/^ root of unity, and the sum 
of all such roots is one of —1,0,1. The map Z/n "Zi/dm induces a map 
(Z/n)* — >■ (Z/dm)*, which is onto, with multiplicity <p{n) / <p{dm) ■ Hence 

(27) lE-^(-r)l< ^^"^ 

k 



Combining Equations 25, 26 and 27 we get 

^^^^ Wp) ^ Wg) ^ MO ^ ^' 

There is no positive integer k such that 0(fe) = 3, and 4'{k) = 2 only for 
k G {2, 3,4, 6}. We conclude that all three roots are 12*^ roots of unity. □ 

Remark. For each choice of A'^, Lemma 3.1 leaves only finitely many cases 
to check. One can finish the proof, for small choices of N, by a quick computer 
search. We did this for N = 8, the case of interest. The argument following 
this remark finishes the proof, in all cases, in a noncomputational manner. 



Lemma 3.2. None ofp,q,r is relatively prime to n. 



Proof. Suppose that p is relatively prime to n. All our congruences will be 
taken mod n. Applying a Galois automorphism of Q[a-'„] we arrange that p = 1. 
Assume first that n is even. Since 5ft(r) < —1 we must have ?R.{uj^+uJn) < and 
^(uJ^ + LUn) < 0. This is only possible if q,r = n/2. But then it is impossible to 
getp+q+r = 0. Similarly, if n is odd then g, r = (nib l)/2. Since p+q+r = 0, 
and n > 2 we must have q = r = {n — l)/2. Here = a;'' is a primitive n}^ 
root of unity. We apply a second Galois automorphism to get q = r = 1. Thus 
?R.{uj'^ + uj^) > 0, a contradiction. □ 

Lemma 3.2 says, in particular, that n must have more than one prime 
factor. Thus, n G {6,12}. Suppose n = 12. Since u;^^a;^a;^ = 1, we see 
that dr divides the least common multiple n' of dp and dq. Lemma 3.2 says 
dp, dq, dr < 12. Suppose n' < 12. Since dp, dq, dr all divide n' we can replace n 
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by n', contradicting the minimality of n. Thus, n' = 12. Hence, {dp,dq) must 
be one of (3,4), (4,3), (4,6) or (6,4). Prom this it is easily seen that dr = 12, 
a contradiction. The case n = 6 is similar. 

Remark. The reader interested only in the parabolic case can skip to the 
next chapter at this point. 



3.4. Near monotonocity of the eigenvector. Let Eg be the C-circle stabi- 
lized by the element gg, defined in Section 3.1. Let E* be the affinely normal- 
ized vector which is polar to Eg. Note that E* is an eigenvector corresponding 
to the norm-one eigenvector As of gg. Since the element r conjugates gg to g~^, 
we have E* = (e(s), e(s), 1). 

Say that a function / : J ^ R is e-monotone if there is a monotone 
function / : J — ^ R such that ||/ — /||oo < £■ 

Lemma 3.3 (near monotonicity) . The functions arge(s) and \e{s)\ are 
-monotone decreasing for s E J. 

Proof. We will show below that |e'(s)| < 3 for all s G J. We explicitly 
evaluate the two functions above for 120001 maximally and evenly spaced 
points in J, and observe that the two functions are monotone decreasing on 
these points, up to a roundoff error of 10~^°. See Section 12 for details of the 
calculation. Using the bound on the derivative, together with the fact that J 
has length less than .54, we get the result of this lemma. 

We turn now to our bound on |e'(s)|. Let Ai and A2 be the functions 
defined in Section 3.2. Let 

X{s) = A2{s)e{s) + ^ 
e(s) 

We first establish the following implication: 

(29) |e(s)| < 1 & \X{s)\ > ^ ^ |e'(s)| < 3 & \X'{s)\ < 20. 

Assume that |X| > 4/3 and |e| < 1 at some parameter s £ J, which we 
suppress from our notation. Note that 

g{e, e, 1) = (Aqc, Aqc, Aq) = (-e, -Aie + A2, -A^e - A{). 

Equating the first and third coefficients respectively gives 

30 Ao = — e = . 

e A-i 

Combining these equations we obtain A^e^ -\- A\e — e = 0. Differentiating with 
respect to s we get 

[^se^ + Xe] + (2^26 + Ai)e' - e' = 0. 
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Basic calculus shows that |^2('S)| < -5 for s G J. If |e| < 1 then the 

bracketed expression above has norm less than 1. Therefore 

(31) 1(2^26 + Ai)e' - e'l < 1 

A2e+- = X. 
e 

\e'{s)\\X-u{,)\<l- u{s) = e'{s)/^{s). 

Since n(s) is a unit complex number and \X\ > 4/3 we see that |X — ii| > 1/3. 
Hence |e'| < 3. Note that |^2(s)| < 3 for s £ J. Hence 

\X'\ < 2\A'2\\e\ + 2\A2\\e'\ + l^l < 1 + 18 + .5 < 20. 

We now know that Equation 29 is true. 

Say that s G J is a good parameter if the two bounds on the left hand 
side of Equation 29 hold at s. Otherwise say that s is a bad parameter. Let 
^0; • • • ) S55 be the set of 56 maximally and evenly spaced points in the interval J. 
Note that J has length ^/125/3 — \/35 < .55. Thus, every point of J is within 
.005 of some point sj. As explicit calculation shows that 

> 1.49; \e{sj)\<.9; j = 0, . . . , 10. 

(See §12 for details of the calculation.) Suppose s G J is some bad parameter. 
Let Sj be the point in our list nearest to s. We can assume that all parameters 
in [sj,s] are good. This interval has length at most .005. Equation 29, says 
that |X| and |e| can each change by at most .1 on this interval. Hence, |^(s)| > 
1.49 - .1 > 1.3333... and |e(s)| < .8 + .1 < 1. This is a contradiction. Hence 
every parameter is good. Our bound follows immediately. □ 

4. Images of the Clifford torus 

4.1. A discreteness criterion. We begin with a discreteness criterion. This 
discreteness criterion is a close variant of the ping-pong lemma and the Klein 
combination theorem. For original sources, see [K] and [Mas]. A formulation 
which is similar to the one here appears in [Mac]. 

Recall that : F — > PU(2, 1) is our basic representation. We have also 
set = Ps(r). Let Hg C be the group generated by h^^g and /i-.s. We say 
that a compressing pair for ps is a pair of subsets {Ui, U2), having the following 
properties: First, d{Ui) = 112- Second, there is a proper subset V C Ui such 
that h{U2) C V for all h G H^. 



Observe that 

2^26 + ^1 = A2e + — = 

e 

Rearranging Equation 31 we arrive at 
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Lemma 4.1. If Ps has a compressing pair then p is a discrete embedding. 

Proof. Any element g ET has the property that 

Ps{g) = hnodo hn-i odo ...ohi 

where hj G Hg. The elements /t„ and hi can be trivial, but all the other hj 
are nontrivial. We perform induction on n. The inductive hypothesis is the 
following: For j = 1,2, 

1. if hn is trivial, then p{g){Uj) C d{V) C C/2; 

2. if hn is nontrivial, then p{g){Uj) C hn{U2) C V. 

When n = 1, there is nothing to prove. 

Suppose the induction hypothesis holds for n — 1. If /in = 1, then 

p{g)Uj = p{d)p{hn-i od...do hi)Uj C d{V). 

If hn I, then 

p{g)Uj = p{hn)p{d o hn-i od...do hi)Uj C p{hn)U2 C V 

as desired. 

As in the ping pong lemma, we have shown that Ps{g) does something 
uniformly nontrivial when g is not the identity element. Hence, ps is a discrete 
embedding. □ 

4.2. The projection. Recall that T is the Clifford torus. — T consists 
of two open solid tori: 

Ui = {{z,w)\ \z\ < \w\}nS^; U2 = {iz,w)\ \z\ > \w\}nS^. 

We will reprove the discreteness result in [GP] by showing that {Ui,U2) is a 
compressing pair for ps as long as s G [0,s]. In this section we reduce this to 
a problem about disks in the hyperbolic plane. 

The condition d{Ui) = U2 is true since d{z,w) = {w,z). Let W C CH^ 
be the closed disk {z = 0} n (CH^ u S^). Let Uw : ^ W he the projection 
Ilw{z,w) = {0,w). For j = 1,2, let Wj = Ilw{Uj). By definition, W2 is the 
open disk of (Euclidean) radius V2/2 centered at (0, 0), and Wi is the annular 
region {(0,u;)| \w\ G (V^/2, 1]}. 




Figure 4.2. 
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Lemma 4.2. (Ui, U2) is a compressing pair for ps if and only if h{W2) 
C Wi for all h G Hg. 

Proof. Using the formula for given in Section 2.3 we see that h±oIlw = 
Uw o h± . Hence, 

(32) hoUw = Ilwoh; h e Hs- 

From this we see that h{U2) C Ui if and only if h{W2) C Wi. If h{W2) t Wx 
then {Ui,U2) is certainly not a compressing pair. 

Suppose h(W2) C Wi for all h e H. Note that h{W2) is a disk and Wi 
is an annulus. Hence, every finite union of the form U^te//' ^(^2) is a proper 
subset of Wi. Here H' is a finite subset of H. The element /i±|i4/ is an isometric 
rotation about a point p± G W, as shown in Figure 4.2. Hence, the product 
h+h- is loxodromic, regardless of parameter. 

The Euclidean diameter of the set h{Ui) tends to as the word length of 
h tends to 00. This fact, together with the finite union property, implies that 
X = U/ieH^(^2) is a proper subset of Wi. Thus, V = U:^;^^{X) IS a proper 
subset of Ui such that h{U2) C V ioi all h e H . □ 

4.3. Patterns of disks. The disks and /i±(VF2) are tangent at the point 
p±. Thus, the disks hh±{W2) and h{W2) are tangent at the point h{p±). The 
points of tangency all lie on the geodesic A, through the points p±. This 
geodesic is the translation axis for hj^h-. Figure 4.3 shows two pictures of the 
situation. 




Figure 4.3. 



550 



RICHARD EVAN SCHWARTZ 



In the second half of Figure 4.3, the complex slice W has been identified with 
the upper half plane model and the geodesic A has been identified with the 
positive imaginary axis. 

If h has odd word length then h is an isometric rotation about a point 
not in W2. Hence, h{W2) C Wi for words of odd length. If h{W2) n 1^2 = 
for some h e H having length 2k then h'{W2) n W2 = for all words h' 
having longer (even) length. The point is that the hyperbolic distance from 
the hyperbolic center of W2 to the hyperbolic center of h'{W2) exceeds the 
hyperbolic distance from the hyperbolic center of W2 to the hyperbolic center 
of h{W2)- Our conclusion: If h{W2) n W2 = 0, for h having word length 2, 
then (Ui, U2) is a compressing pair at the parameter in question. 

By symmetry, h+^sh-,s{W2)nW2 = if and only if h^^,,h+^s{W2)nW2 = 0. 
A calculation, which we omit, shows that the distance between the hyperbolic 
center of W2 and the hyperbolic center of /i+^2^-,s(1^2) is monotone decreasing 
in s, and that the two disks are tangent when s = s. This establishes the 
discreteness result in [GP]. 

Remark. Our discreteness criterion parallels [GP, Lemma 3.2]. Our Lemma 
4.2 parallels [GP, Lemma 3.3]. The analysis in this section parallels the auxil- 
iary lemma used to prove [GP, Lemma 3.3]. The omitted calculation parallels 
the calculation in [GP] concluding their discreteness proof. Our Clifford torus 
is contained in the ideal boundary of the Dirichlet polyhedron considered in 
[GP]. The map Hw corresponds to the map He used in [GP]. 

4.4. Another point of view. The reader may think, at this point, that our 
proof above fails by accident. Let us give another point of view which shows 
that the proof must fail for structural reasons. 

The fixed point pj of gj contained in the Clifford torus, by symmetry. 
The point here is that the element r preserves the Clifford torus and conjugates 
Qs to its inverse. Let Q be the complex line tangent to S"^ at pj. The element gj 
stabilizes Q, and rotates this plane nontrivially. Q is transverse to the Clifford 
torus T. This means that gs{T) must intersect T transversely. For s sufficiently 
close to s the same phenomenon must occur, by continuity. One may think of 
s as the cutoff for this phenomenon. When s € (s, s], the element Qs does not 
have the translational strength to compensate for its rotational component in 
a direction transverse to T. In short, gg does not have the strength to pull T 
off of itself. 

5. The hybrid cone 

5.1. The basic construction. Suppose that (E,p) is a pair, where £^ is a 
C-circle and p e E is & point. 
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Lemma 5.1. Suppose x G — E. There is a unique H-circle 7 = 
^{E,p; x) such that x E j and p € 7 and ^r\{E — p) 

Proof. We normalize by a Heisenberg stcrcographic projection so that 
E = Eq, the center of TC, and p = 00. In this case, there is a unique level 
Heisenberg R-circle containing x. This is 'y{EQ,oo]x). □ 



Q(£o,;7;x) 



Figure 5.1. 

Let x) be the arc of ^{E,p; x) which connects X to p and which is 

disjoint from E — p. li S C — E we define 

(33) n{E,p;S)= \JniE,p;x). 

xes 

We call Q{E,p;S) the hybrid cone of S with respect to {E,p). We call the 
individual R-arcs, in the union above, the foliating arcs. 

Our construction is natural: n{T{E),T{p);T{S)) = T{n{E,p; S)), when 
T is the restriction to of an isometry of CH^. 

We define the endpoint map 

(34) fn{x)=-1{E,p■x)r^{E-p). 

Remark. When fa{S) is a single point, ^7 is a subset of a spinal sphere. 
(See [G, Ch. 5] for the theory of spinal spheres.) Thus the hybrid cone con- 
struction generalizes the spinal sphere construction. 

5.2. Hybrid disks and hybrid sectors. Let {E,p) be as above. If S is 
a C-arc, let S be the C-circle which contains S. A hybrid disk is a hybrid 
cone of the form il{E,p;S). A hybrid sector is a hybrid cone of the form 
Q,{E,p;S), where 5 is a proper nontrivial C-arc. We will sometimes write 
n{E,p; S) = n{E,p; S). 

Lemma 5.2. Let = il{E,p;S) be a hybrid sector. Q, — dO, is an ana- 
lytically embedded open solid triangle, ft — p — C is an analytically embedded 
punctured open disk. 
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Figure 5.2. 

Proof. We will consider the case of hybrid sectors, the case of hybrid 
disks being similar. We normalize by a Heisenberg stereographic projection so 
that E = Eq, the center of H, and p = oo. In this setting 0(£'o,oo;5) — oo 
is a surface with boundary, ruled by rays. The ruling rays are arcs of level 
Heisenberg R-circles. The endpoints of the rays are contained in S. The 
linking condition implies that the rays are pairwise disjoint. The rays certainly 
vary analytically with points on S. Our lemma is clear from this description. 

□ 

When i}{E,p; S) is normalized as in the previous lemma, we say that it is 
in standard position. Figure 5.2 shows a picture of a standard position hybrid 
sector. The point oo is not shown. If continued in the other direction, all the 
foliating arcs of a standard position hybrid sector intersect the center Eq. This 
accounts for the radiating appearance of the top view. 

5.3. Canonical parametrization. In this section we describe a parametriza- 
tion for hybrid sectors. The formulas for the parametrization are only impor- 
tant for the proof in the sense that our code actually uses them. 

Let {E, p; S) be as above. Let q be the point of E diametrically opposite p. 
(This makes sense because is a round circle.) Let C, ^ he some auxiliary 
point chosen so that p, q, C lie on a common R-circle. Let {p, q, ^} be an 
equalized choice of lifts. Let E* be the polar vector to Eg, chosen so that 
{C,E*) = {p,q). 

Let 5 be a C-arc as above. Let Si and ^2 be the endpoints of S, and let 
S3 be the midpoint of S — S. Referring to Equation 10 we define 

(35) St = Css{Si,S2;t), St = j^:^^^, te[0,l]. 
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St is the lift of St such that {St,p) = 1. Now define 

(36) At = ^{{q, St)) p + iq, At = 6(1*). 

Note that At lies on the complex slice determined by p and q. It is easily 
checked that {At, At) = and A{p,St,At) = 0, so that At = fn{St). With 
reference to Equation 13, the following map is a paranictrization of ^{E,p; S): 

(37) n{E,p;S;t,u) = RA,ip,St;uy, {t,u) e[0,lf. 



5.4. The elevation map: Basics. With notation as above, the formulas in 
this section are only relevant to the proof in that we actually use them in our 
code. 

Recall that it{z,w) = (z,S(w)). Define 

(38) Y^ = 7roeoY^; Y^{x) = {{x, E*), {x,q), {x,p)). 

It is easy to see that the vector V{z,w) = wp + q + wE* is null if and only 
if 3ft(w;) = \z\^. Evidently G o Y{V{z,w)) = {z,w). Thus 6 o y o 6"^ maps 
S^ — p into the boundary of the Siegel domain. In short, is a Heisenberg 
stereographic projection mapping n{E,p;S) to standard position, normalized 
so that Y^{p) = oo and Y^{q) = (0,0) and y^(C) = (1,0). 

Let S^ be the unit circle in C. Let T = x R. We define 

(39) ^{z,t) = {z/\z\,ty, 7rn{z,t) = t. 

Next, the (^-elevation map ^ : S^ — {p} — s- T is defined as follows: 

CoY(;{x) ifxeS^-E 

S^ X TT-R, o Y^{x) iixeE-p. 

Technically, is a map from S^ — p to subsets of T. This detail is absorbed 
into the statements of our estimates. 

Lemma 5.3. ^{niE,p; S) - {p}) = ^(5). 



Proof. Let a: G be a point. By construction, Y (yi(E , p; x)) is a ray of a 
level R-circle. Thus, ^ is constant on Y {^1{E , p; x)) . □ 

Say that a barrier for {E,p) is a pair (/, ^), where / : S"""^ ^ R is a 
continuous function. Let Aj be the graph of / in T. Let Aij and A2J be the 
two components of T — Aj. We will take A2J to be the component above Af, 
as shown in Figure 5.3. In this figure, 5^ x R is mapped into the plane via the 
map {z,t) — ^ {aig z,t). 
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Figure 5.3. 



Lemma 5.3 immediately implies the following result: 

Lemma 5.4. Let Q = Q{E,p;S) be a hybrid sector. Suppose that 7 G 
PU(2, 1) stabilizes {E,p). Suppose C Aij and ^'(7(0)) C A2,/. Then 

nnjin) = {p}. 

Remark. At this point, the reader interested in the parabolic case of the 
Main Theorem can skip to the next chapter. 

5.5. The elevation map: Fine points. 

Lemma 5.5. Suppose that a is an R-arc, contained in an ^-circle a. If 
a intersects E in two points, neither of which is p, and a intersects E in one 
point, then ttr o"^ is monotone on a. 



a 




top view 



E 



Figure 5.4. 
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Proof. We set Y = Y^, and recall that ^' = S,oY, as above. Also, ■Kc°Y{a) 
is a lemniscate centered at the origin and Y{a) is an f -integral lift of this 
lemniscate. The slope of this lift is only at the two points which cover 0. In 
particular, the slope of Y{a) vanishes only at the endpoint. (See [G, p. 141] for 
an explicit parametrization.) Hence, the height function is monotone on ^(q;). 

□ 

Corollary 5.6. Let Q{E,p;S) be a hybrid disk. Suppose 7 G PU(2, 1) 
is such that "y{E) = E but ^{p) 7^ p. Suppose that ^ o ^{dVi) C A2 and 
ttr, 0^*0 7(p) > max/. Then ^{a) C A2J for some foliating arc a of 'y{il). 

Proof. The curve ^ o 'y{dQ,) is a simple closed curve on T which, by 

hypothesis, is contained in A2J. This curve is a generator for the homology of 
T. From this it follows that there is some x G jidi^) such that ttr o vJ/^x) > 
max/. Let a be the R-arc of 7(^^) which connects x to p. By Lemma 5.5, we 
have TTn°^{y) > max/ for all y E a. This obviously implies that '^{a) G A2,/. 

□ 

5.6. Tame hybrid disks. Say that a hybrid disk Cl{E,p; S) is tame if the 
surface ^ — dfl — p is never tangent to the distribution 6. 

Lemma 5.7. Let n{E,p;S) be a hybrid disk. Let E* and S* be polar 
vectors to E and S. Ifd{E*,S*) < 1/4 then Q,{E,p;S) is tame. 

Proof. Modulo PU(2, l)-equivalence, there is a one-parameter family of 
inequivalent hybrid disks. To get a complete family, we normalize so that 

E* = {0,1,0); 5* = (l/a,l,l); p=(-l,0). 

Here a G [0, 1). We compute 6{E*, S*) = a^. The map Bq, defined in Equation 
7, maps to standard position. The points 

/ N f 2a a2-l\ 

Pi = (0,1); P2= { . 2 ' 1 _L 2 
\ 1 -I- a"^ 1 + a"^ J 

belong to S. By symmetry, Bo(-S') intersects R x {0} at Bo(pi) = (1,0) and 
Bo(P2) = (a— 1,0) and has center of mass c = (a,0). The point on 7rcoBo(n) 
having minimum norm is ttc o Bo(j'2) = a — 1. If a G [0, 1/2) then every point 
of TTc o Bo(ri) has norm greater than a. 

Let Tx and E^ be the tangent plane and £^-plane at x G Bo(fi — iS) — 00. 
Given any plane H G we define the slope s{lV) = sup^gn ^R(^)/Kc(^'')|• 
From Equation 8 wc have s(Ex) = 2|7ri(x)| > 2\a\. Now, Bo(ri) is ruled by 
horizontal rays, which have endpoints on S, and which lie on lines that intersect 
{0} X R. This geometry implies that s{Tx) < s{Ec) = 2\a\ < s(Ex). Hence, 
Sx 7^ Ex. □ 



556 



RICHARD EVAN SCHWARTZ 



a-\ 




R 



Figure 5.7. 



5.7. A disjointness criterion. 

Lemma 5.8. Suppose that I C R"" is some connected open set. Suppose, 
for t E I, that Q{E,pt;Ct) is a smoothly varying family of tam,e hybrid disks 
and ^l{E,qt, St) is a smoothly varying family of hybrid sectors or hybrid disks. 
Suppose that 

1. For each t £ I, pt qt. 

2. Ct n ^{E, qt] St) = for all t G /. 

3. St n ^{E,pt; Ct) = for all t G /. 

4. For some parameter s E I, one of the foliating arcs of Q,{E,qs, Ss) is 
disjoint from il{E,ps, Cg). 

Then n{E,pt, Ct) n Q,{E, qt, St) = for all t G /. 

Proof. Let N denote the union of all foliating arcs of Q,{E,qt,St). The 
parameter t ranges throughout / in this definition. We equip with the 
Hausdorff topology: Two elements of iV are close if they are contained in small 
tubular neighborhoods of each other. Obviously, is connected. 

Let N' C N denote those arcs /? such that /? is a foliating arc for Q{E, qf, St) 
and f3 is disjoint from n{E,pt, Ct). To establish this lemma, it suffices to prove 
that N' = N. Note that N' is nonempty, due to statement 4 above. Clearly, 
A^' is open. 

To finish the proof we show that N' is closed. Let /32... be a sequence 
of elements of N' which converge to some (3 e N. Let t be the parameter 
associated to /?. 

1. Statement 1 above implies that /Jflpt = 0. 

2. Statement 2 implies that /? n Ct = 0. 

3. Statement 3 implies that dp D n{E,pt, Ct) = 0. 
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The only possibility is that (3 — dfi intersects Q,{E,pi,Ct) — Ct — pt- If this 
intersection is transverse, then nearby arcs, including some of the /?„, would 
intersect nearby hybrid disks. Hence (3 — 8(3 is tangent to fl{E,pt;Ct) — Ct — pt 
at some point x. 

There is a foliating arc of fl{E,pt,Ct) which contains x. Call it a. Since 
n{E,pt,Ct) is tame, the plane E^ of S intersects the tangent plane to 
Q,{E,pt; Ct) transversely, in a line L. Both a and (3 are tangent to Q,{E,pt, Ct) 
at X, and both a and (3 are £^-integral. Hence, a and (3 are both tangent to L. 
The important conclusion here is that a and (3 are therefore tangent to each 
other. 

We normalize so that Cl{E,pt;Ct) is in standard position. Also, Trc{a) 
is a straight line segment, contained on a ray through the origin and /? is a 
Heisenberg R-arc, contained in a Heisenberg R-circle p. The curves 7rc(Q!) 
and nc{(3) are tangent at ircix) 7^ 0. 

There are two cases to consider. Suppose first that 7rc(/?) is a straight 
line. In this case, 7rc(a) and 7rc(/3) are the same line through the origin. This 
implies that a = (3. Hence da fl /3 7^ or da n /3 7^ 0. Either possibility 
contradicts our assumptions. 

We conclude that TTciP) is a lemniscate. Since intersects Eq twice, 
TTciP) is centered at the origin. However, a straight line (7rc(Q;)) through the 
origin cannot be tangent to a lemniscate centered at the origin (7rc(/3)) at a 
point (7rc(a;)) other than the origin. This is a contradiction. □ 



6. Dented tori 

6.1. Technical calculation. 

Remark. For the reader only interested in the parabolic case, the first two 
sections of this chapter can be replaced by a single calculation at the parameter 
s. One can just skim this material. 

Recall that E* = (e(s), e(s), 1) is the eigenvector to gg, as in Section 3.4. 
Let S be the circle of radius \/2/2. Let pi{s) and qi{s) be the points such that 
the tangent lines to S at pi{s) and qi{s) contain e{s). That is: 



(40) Pi{e) = ^— , P2{e) 2e 

Figure 6.1 shows a picture which is visually accurate for all s G J. 
The following result implies that pi{s) never exits the shaded region of 
the picture and qi{s) never enters the shaded region. 

Lemma 6.1. For all s e J, argpi(s) G [7r/12,7r] and argp2{s) G 
[7r,27r + 7r/12]. 



558 



RICHARD EVAN SCHWARTZ 




Figure 6.1 

Proof. Define 9i = arg (e(s)) + 10-^ and 62 = arg(e(s)) - IQ-^. Also define 
fi = |e(s)| + 10~^ and r2 = |e(s)| — 10~^. Finally, define eij = riexp(^j). 
From the near monotonicity lemma, and basic geometry, pi lies in the interval 
bounded by the outermost points Pi(ejj) and qi lies in the interval bounded 
by the outermost points qi{eij). Explicitly calculating, we obtain this lemma. 
Indeed, the points pi{eij) and qi{eij) are positioned almost exactly as those in 
Figure 6.1, independent of indices. □ 



6.2. The axis lemma. Define 
(41) a = tan(7r/2 - 7r/12) =2 + ^3. 

Let Act be the corresponding Clifford triangle, on the Clifford torus T. Let 
dhA be the arc of the horizontal C-circle which bounds A. Let d^A be the arc 
of the vertical C-circle which bounds A. Let d^A be the arc of the diagonal 
C-circle which bounds A. Thus, we have dA = d^A U d^A U d^A. 



Lemma 6.2 (axis lemma). Eg CiT consists of two points. One of them, 
Ps, is contained in Aa-. The other point, qs is not. Also Es links each of the 
three C- circles in dA^^. 



Proof. We repeat the notation from Section 6.1. Let S be the disk of 
radius y/2/2 centered at the origin of C. Let p\{s) and q\{s) be the two points 
of dS such that the tangent lines through these points contain e(s). Define 

(42) Ps = {Pi{s),Pi{s)); qs = {qi{s),qi{s)). 

We claim that E's H T = pgU qg. By consruction, Ps,qs £ T. Recall that 
R = {z = w} n 5'^. Note that Ps,qs S R. Since i? is a round R-circle, it 
is contained in the real slice R = {z,z} C C^. The map — > C given by 
{z,z) — z carries R to the disk S and carries geodesies in R to straight lines 
in C. Hence, the complex lines tangent to at Ps and qg both contain E*. 
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This is enough to conclude that E* is Hermitian perpendicular to Ps and Qs- 
Hence, EsC\R = Ps^Qs- If some other point Og lies in EgClT, then so does r{os) 
which is contained in the same diagonal chain as Og- This would force Eg to 
be a diagonal C-circle, which it is not. The result of Section 6.1 is equivalent 
to the statement that Ps & A and gr^ A for all s G J. 
Calculations show that 

(43) 5{E^, ddA*) = .4; S{E^, dhA*) = 5{E^, d^A*) = .042938 . . . . 

Here ddA* is a polar vector to the C-circle containing ddA. Likewise for the 
other vectors. 

Let X be any one of d,v,h. The calculation above shows that Ej links 
the C-circle containing d^A. Since Eg never intersects the C-circle containing 
dxA for s G J, we see that Eg always links this C-circle. □ 

6.3. Main construction. We define 

(44) n^is) =n{Eg,pg; d^A); x = h,v, d. 

For any subset S C let dS = d{S). Here d{z,w) = {w,z). Next, define 

(45) J^a(s)= IJ n^is); n{s) = nA{s)UdnA{s). 

x=h,v,d 

By construction, il{s) is the union of 6 hybrid sectors. Next, we define 

(46) S = T - A - dA; Z(s) = n{s) U S. 

We call Z(s) the dented torus. Z(s) has been constructed by cutting out two 
Clifford triangles from T and gluing back in 6 hybrid sectors. Figure 6.3 shows 
some of these sets schematically. 




Figure 6.3 
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d preserves Z(s) by construction, r preserves A, dA, E, Es, Ps, dEg and 
dps, all the sets which are relevant to the construction. Thus Z(s) has 4-fold 
dihedral symmetry. In particular, the points ps and Qs, which are fixed by r, are 
diametrically opposed on the circle Eg. Hence, the canonical parametrization 
of ^x{s), for X = h,v, d, uses the point qg as the extra point. 



Since PsjQsiC £ the ^-elevation map is well defined for the pair 
{Es,pg). By construction ^'5(^^2,(5)) is in standard position, for x = h,v,d. 



7.1. Five disks. Let W be the complex slice used in Section 4. Let 
W* be the complex line {z = 0} which contains W as its unit disk. Let 

= {{0,w) G W*\ ^{w) > 0}. Define Wl similarly. Let C be the 
translation axis of the element hj^^h^^. 

We define five disks which play a key role in our proof. 

1. Let Wa{s) C W* be the disk such that Ag c dWA{s). 

2. Let W+{s) = h+-s{W2). 

3. Let Ws{s) = h--g{W'i). 

4. Let W^{s) = h+;g{W+*). 

5. Let W^{s) = h+;g{W-*). 

Let Wabc{'s) be the union of these five disks. Figure 7.1 is a fairly accurate 
picture (compare Figures 13.1-13.6). Henceforth we suppress the parameter s. 



Define 



(47) 




7. The discreteness proof: Parabolic case 




Figure 7.1 
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Let H2 be the set of words in H which have length 2. 

Lemma 7.1 (the five disk lemma). For all nontrivial h & H — H2, h{W — 

Wabc) C Wabc and h{W2) C Wabc- 

Proof. We first deal with W — Wabc- Suppose first h has odd length. By 
construction, h preserves the axis A, and interchanges the two components of 
W - A. Therefore, 



By construction, W — Wq — Wq is a fundamental domain for the action 
of h^h-h-\.h-. Thus, if h has word length at least 4 we have 

(49) h{W - Wabc) C h{W -W^ - W^) c U VF^ c Wabc- 

We now turn to the images of W2- An explicit calculation shows that 

W2 n W^ = 0. Thus, W2 CW -W^ CW - Wabc- Equation 49 now imphes 
that h{W2) C Wabc when h has word length at least 4. If h has word length 
1 then h{W2) C Wb, by definition. 

Suppose h has word length 3. By symmetry we can assume that h = 
hUi^h^. Since n = and W2 n = 0, we have W2 - Wa CW- 

Wabc- Thus, we are reduced to showing that /i+/i_/i+(W2nM^_A) C Wabc- Let 
X = h+{W2 n Wa)- Let X± = X n W^. By construction, h+h-{X+) C W^. 
The region X_ if nonempty, is contained in Wa, by convexity, as shown in 
Figure 7.2. Thus, 



(48) 



h{W - Wabc) c h{W - Wa) cWaC Wabc- 



h+h^{X^) C h+h.{WA) =WaC Wabc- 



Ah in ah, /i_,_(W2 n W a) = h+h-{X) c Wabc, as desired. 



□ 




X_ inhere 



Figure 7.2 
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7.2. All words but two. Let Z = Z(s) be the dented torus constructed in 
the previous chapter. Recall that Z = U S. 

Lemma 7.2 (estimate 1: parabolic case). H Wabc = 0- 

Proof. See Section 12. □ 

Corollary 7.3. Ilw{^) n Wabc = 0- 

Proof. Let r and d be as in Equation 18. This corollary follows from the 
facts that J7 = U d^lA, rd commutes with XI^^, rd{WABc) = Wabc and 
rd{Vti^) = d^A- □ 

Lemma 7.4. If h e H - H2 is nontrivial then h{Z) n Z = h{E) n S. 

Proof. By the five disk lemma, and Corollary 7.3, we have h{Il\Y{Q)) C 
Wabc- Hence nH^(il) n Uw{h{Q)) = 0. Pulling back, we have n n = 0. 
By the five disk lemma, we have h{W2) C Wabc- Hence, by Estimate 1, 
we have Ilw{^) n h{W2) = 0. Pulhng back by n,7y^, and using the fact that 
nvK(2) C W2, we get O n h{E) = 0. Since this is true for all nontrivial words 
he H - H2we have n n h~^{E) = 0. That is, h{n) n S = 0. This result now 
follows from the fact that Z = ft U E. □ 

Let X be the attracting fixed point of h^h^. This point is one of the two 
intersection points of the axis A with dW. In particular, 

X = Uw{x) c Wa c Wabc- 

By Estimate 1, we have x Z. Let Ui he the component of — Z which 
contains x. 

Lemma 7.5. Ifye Wabc - W2 then U^{y) c Ui. 

Proof. Note that Wabc ~ W2 is path connected. Hence, there is a path 
5i C Wabc ~ W2 which joins y to x, the fixed point of hj^h^. We can find 
a path 82 C S'^ such that Iiw{^2) = ^1- One endpoint of 82 is necessarily x, 
since Jl'^{x) = {x}. We can make the other endpoint any point of Ii'^{y) we 
choose, since the fibers of Hw are connected. By construction, 82 d — Z. 

□ 

Corollary 7.6. Both fixed points of h+h- are contained in Ui. If a is 
a fixed point of dh+h-d and h E H then h{a) G Ui. 

Proof. For the first claim, the fixed points of h+h- project to 9VFn W^i C 
Wabc ~ W2- For the second claim, the fixed points of h^h^ have zero first 
coordinate. Hence, a has zero second coordinate: Hpi/(c7) = (0,0) Wa- If h 
has odd length then Ilw{h{a)) = h{0, 0) C Wa — W2- If h has even length, at 
least 4, then Ilw{h{a)) = /i(0, 0) C Wq, which is disjoint from W2. If h has 
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length 2 then /i(0, 0) is contained in dW^ . We have already seen that W2 is 
disjoint from this disk. So, in all cases, ^(0, 0) G Wabc — W^2- Pulling back by 
liw, and applying Lemma 7.5 finish the proof. □ 

Lemma 7.7. C/i ^1/2. 

Proof. Let y be the fixed point of the parabolic element h+d. Recall that 
iJ+ is the horizontal chain on the Clifford torus which contains the horizontal 
arc bounding A, the Clifford triangle used in the construction of Z. By con- 
struction y G Hence, Ilw{y) is the fixed point p-^. G W2 n Wa of h+. In 
particular, Uwiu) C. Wabc- 

Let X be the fixed point of h^h-. There is a path Sq, connecting x to y, 
such that So — {y} G Ui, and So is transverse to S at y. The path Si = SoUd{do) 
joins X E Ui to d{x) G U2- By construction, Si intersects H transversely in a 
single point and is disjoint from Z — 3 = Q. If Ui = U2 then there is some 
path 6[ which connects x to d{x) but which avoids Z. The path 7 = 5i U 6[ is 
a closed loop in which intersects E once, transversely, and is disjoint from 
Jl. Hence 7 and dO, = dE are algebraically linked. 

Homological considerations say that 7 links dA or 7 links ddA, or both. 
Suppose, without loss of generality, that 7 links dA. Since is the continuous 
image of a disk, attached to dA by a degree one map, we see that 7 and Qa 
have nontrivial intersection. This is a contradiction. □ 

Lemma 7.8. Let h G H — H2 be nontrivial. The component Ui contains 
all but one component of — /i(Z). 

Proof. By Lemmas 7.3 and 7.4, the intersection X = /i(Z)nZ is either the 
empty set, or one of the two C-circles H±. In either case, X neither disconnects 
Z nor h{Z). Basic topology now implies that one component of — Z, which 
we call the special component, contains all but one component of — h{Zi) 
(and vice versa). 

We now show that Ui has nontrivial intersection with more than one 
component of — h{Z), meaning that Ui is the special component. Since h 

interchanges the fixed points of h+h-, both of which arc contained in Ui, we 
sec that h{Ui) n ?7i / 0. From Corollary 7.6 we see that Ui also contains h{a). 
Thus h{U2) n [/2 / 0. □ 

Lemma 7.9. For all nontrivial he H - H2, h{U2) C Ui. 

Proof. We already know that Ui contains all components of S*'^ — /i(Z), 
except for one. We just have to rule out the possibility that h{U2) is this one 
exceptional component of — h{7i). Assuming to the contrary, we get h{Ui) C 
Ui. The element rd conjugates h to h~^ and preserves Z. Furthermore, since rd 
interchanges the fixed points of h+h-, we see that rd{Ui) = U\. If h{Ui) C U\ 
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then 

h-^{Ui) = {rd)h{rd){Ui) = {rd)h{Ui) C rd{Ui) = Ui. 

This impUes that h~^{Ui) n f72 = 0. By Corollary 7.6 we have h'^{Ui) n C/2 D 
{a} 7^ 0, a contradiction. □ 

Remark. Computer evidence suggests that Z(s) is a tamely embedded 
topological torus and, in particular, that — Z(s) consists of exactly two 
connected components. If this is true, our argument above could be somewhat 
simplified. While it is possible to convert the computer evidence into a proof, 
it is quite difficult to do so. We settle for the slightly complicated arguments 
above to avoid the extremely complicated arguments needed to analyze the 
topology of Z(s). 

7.3. Words of length 2. Recall that T = x and : - p ^ T is 
the ( = (l/\/2, l/\/2) elevation map associated to {E,p). Define the barrier 
function: 

(50) f{u)=t + .2m{zu)- (z^2^) = * 05(C). 

In defining z, the branch of ^ is chosen so that < arg(z) < arg(2;^). Let 
A and Kj be sets associated to the barrier (/,*)• We found this function 
experimentally. 

Lemma 7.10 (estimate 2: parabolic case). For the parameter's: 

1. ^'(HUdJ^A) C Ai(/). 

2. ^'o5(dJ7AUS) c A2(/). 

Proof. See Section 12. □ 
Corollary 7.11. For h e H2, h{Ui) c Ui. 

Proof. It suffices to consider the element hj^h-, by symmetry. Note that 
/i+/i_(Z) = h+h-d{Z) = g{7,). Prom Estimate 2 and Lemma 5.4 we have 
g{Z) n Z = {p}. Hence /t+/i_(Z) n Z = {p}. Since p disconnects neither Z nor 
/i_l_/i_(Z), the same argument as above shows that /i+/i_(C/2) C Ui. □ 

All in all {Ui, U2) is a compressing pair for p, and p is a discrete embedding. 

8. The year lemma 

Recall that J = [s, s]. We partition J into 365 intervals of equal size. The 
set Jo of midpoints of these intervals is given by 

(51) (1 _ + fc = l,...,365. 
\ I ^ 730 730 

Let 7fe be the partition interval containing s^. 
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The geometry of Z(s) varies slowly with s G J. However, this variation is 
extremely difficult to estimate directly. In this chapter we replace the family 
of sets {Z(s)| s G J} by a family of sets {Z{sk,t)\ t e h; k = I, . . . ,365}. The 
set Z(sfc,t) is a replacement for Z(t), whose geometry is essentially the same 
as the geometry of Z(sfc). In this way, we only have to deal with the geometry 
of the surfaces Z(sfe), for fc = 1, . . . , 365. 

8.1. Compatibility maps. Let s = s^, and t & 1^ he some other param- 
eter. Let {Es,Ps) be the pair used to define Z(s). There are many maps 
(f)st € PU(2, 1) such that (f)st{Es,Ps) = {Et,pt)- In this section, we construct a 
canonical one which varies contim.ioT.isly with the parameters. 

Let i^s be the midpoint of the shorter of the two arcs of R which connects 
Ps to Qs- Likewise define /if Let = 0(i?^) = (e(s), e(s)). Likewise define rt- 
We define (f)st = Pst ° <^st, where ast,Pst £ PU(2, 1) stabilize R and 

1. Ust is a rotation such that a(/Lis) = //t, 

2. Pst is a pure translation along the line which joins ast{rs) to rj. 

By construction, (pst commutes with the map r, which stabilizes R. 

Let p be the round metric on S^. For any map 7 G PU(2, 1), define 
I7I = sup^£53 p(a;, 7(5)). In particular, define Est = {(pstl- As t — > s, we have 
Sst ^ 0. 

For any set X C let be the £ tubulr neighborhood of X, as measured 
in p. Define 

(52) B{s,t) = {dAY; £ = 
B{s,t) is a thin solid torus which surrounds dA. 

8.2. Main construction. Define Q{s,t) = ^A{s,t) U d^Ais,t). Here 

(53) n'^{s,t) = (l)st(.nA(.s)); nA{s,t) = n'^{s,t)UB{s,t). 
Finally, define 

(54) Z{s,t) = Eun{s,t). 
For later use, we define the elevation map 

(55) ^',t = ^',U 

Here is the C-elevation map associated to the {Es,Ps)- By construction 
is the (^'-elevation map for {Et,pt)- Here (' = 0st(C)- 
By construction, 

n'^{s,t) = n{Et,pt;(l)stdA). 
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Thus, "ifst maps the hybrid sectors in Q,'^[s,t) to standard position. In partic- 
ular, the hybrid sectors in fi'^(s, t) are as well adapted to the element gt as are 
the hybrid sectors in the discarded Q^it). 

At the same time Q'^(s,t) is a PU(2, l)-image of ^^is), and thus has es- 
sentially the same geometry. The object Z{s,t) retains the dihedral symmetry 
of Z(s). Finally, the set B{s,t) seals up the gaps created by the discontinuous 
nature of our construction. 

We will not actually deal with Z(s,t) directly. Rather, we will estimate 
the mismatch between Z(s,t) and Z(s). 

Define W^{s,t) = (psti^dis))- The remainder of the chapter is devoted to 
proving the following result: 



Lemma 8.1 (the year lemma). Let a = 
and t G Ik- 

1. Est < a. 

2. gt{EUdnA{s,t)) C {gs{EUdn{s))K 

3. gt{dn'ais,t))c{gs{dnd{s)))'. 



.00046 and b = .00221. Let s = Sk 



To avoid redundancy, we will just prove Statements 1 and 2. Statement 3 
has the same proof as Statement 2. 

8.3. Computational ingredients. 

Lemma 8.2. Let Sk be as above and let t G Ik- There is a constant 
such that 

1. 8||e(sfc)| - |e(t)|| + I arge(sfc) - arge(t)| < ak < a. 

2. 4afe -I- 2\sj -t\<b. 

Proof. Let s = sq, si, . . . , S730 = s be the set of 731 maximally and evenly 
spaced points in Jq. A straightforward calculation, detailed in Section 12, 
shows that, for successive values Sj, Sj+i we have 

1. 8(||e(sj)| - |e(sj+i)|| -|- 10"^) -|- | arge(sj) - arge(sj+i) -|- 10"^! = aj < a. 

2. Aaj + 2\sj — Sj+i| < b. 

Each partition interval has the form [sk-i,Sk+\], with midpoint Sk- Here 
k is odd. The two items above, combined with the near monotonicity lemma, 
give the result. □ 
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Let Hz : ^ C be projection onto the first factor. Referring to Equation 
23 define 

(56) Sit = {x e S^\ sup \Uzix) - Pois)\ > 1/k}, Po{s) = 

seJ ^2(.sj 

The Lipshitz lemma and the variation lemma below explain the signifi- 
cance of this set. 

Lemma 8.3 (part of Estimate 1). For all s G Jo, (S U d^Ais))^ C S3. 
Proof. See Section 9, Estimate 1. 

For later use, let Fj. C T,k x denote those points (x, y) which can be 
joined by a geodesic segment— i.e. an arc of a great circle— which is contained 
entirely in S^. 

8.4. Proof of Statement 1. Statement 1 of Lemma 8.2, together with the 
following lemma, finishes the proof of Statement 1 of the year lemma. 



Lemma 8.4. est < I arg(e(s)) - arg(e(t))| + 8||e(s)| - |e(i)||. 



Proof. Recall that (pst = Pst ° ctst- Here agt is a rotation which maps the 
midpoint /ig to the midpoint //f. Now, /is is contained in the line through the 
origin which also contains the point Q{E*) = (e(s),e(s)). The corresponding 
statement holds for fit- Prom this it follows that 

\ast\ < I arg(e(s)) - arg(e(t))|. 

It remains to show that \/3st\ < 8||e(s)| — |e(i)||. From the near monotonicity 
lemma, and from explicit evaluations at the endpoints of J, 

57 sup — - < 8. 

seJ 2|e(s)|^ - 1 

Let R be the real slice bounded by the round R-circle R. Note that R = 
{{z,z)} n CH^. We put the Euclidean metric on R. By construction, the 
points 

(0,0); Xi = ast{eis),e{s)); X2 = {e{t),e{t)) 

are all contained in a single line L. The element (3st is a translation along L 
which maps pi = ast{e{s) ,e{s)) to p2 = {e{t),e{t)). Note that 

\\pi\\ = V2\e{s)\; IIP2II = V2|e(i)|; - P2II = V2||e(s)| - |e(t)||. 

Combining this information with Lemma 2.3 and Equation 57 gives our esti- 
mate on |/3st|. □ 
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Corollary 8.5. For any point x G H U dQ,A{sk,t), there is some y G 
S U dCl/^{sk) such that p{x, y) < a^. 

Proof. Let s = s^- This follows from the fact that (j)st moves every point 
of less than a^. Hence, dQ,'^{s,t) is contained in the a^-neighborhood of 
drt/^{s). By construction, B{s,t) is contained in the a-tubular neighborhood 
oiEUdnAis). □ 

8.5. Proof of Statement 2. 

Lemma 8.6 (the Lipschitz estimate). For all {x,y) G -Pfe, and all k>0, 
p{gs{x),gs{y)) ^ 2k 



sup 

{x,y)eFk 



< — V27 + 31fc2. 



p{x,y) 

In particular, the quantity on the right-hand side is less than 4 for k = 3. 



Proof. From the formula in Section 3.2 we have 
w Aiz — A2 



g{z,w) 



{gi{z, w),g2{z, w)). 



^A2z + Ai A2Z + A1. 

Note that |^2(s)| > 3^/3/2 for s G J. Let Pq = Pq{s) be as in Equation 56. 
We compute 



\dz9i{s)\ 
1^^52(5)1 
\dw{9i] 



A2W 



{A2Z 


-ArY 


Al- 




{A2Z 


-A^f 




1 



< 



< 



< 



< 



2fc2 



\A2\\z-Pq\^ - 3^3' 



< 



|A2|2|z-Po|2 - 27 



< 



2k 



\A2\\z-Pq\ - 3^3' 
dw92 = 0. 

The bound in this lemma follows from integrating the differential bound 
\\d9{v) 



□ 



< \^\dz9i? + \dw9i? + \dz92? + \d^g2? 

\\v\\ ' 

along the arc of the great circle joining the two relevant points. 



Lemma 8.7 (the variation estimate). For all s,t & J and x G S^, 

Pi9six),9tis)) < y|s-*l- 
In particular, the quantity on the right-hand side is less than 2\s — t\ for k = 3. 
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Proof. Recall that 

g = &ogoe^^; e{zi,Z2,Z3) = {zi/z3,Z2/z3y, &~'^{z,w) = {z,w,l). 

Let X = {z,w) G S^. Let X = Q-^{X) = {z,w,l). 

Inspecting the formula for g, we see that the third coordinate of g(X) has 
norm at least 

q /q 

\Ai + A2z\ > \A2\/k > 

Here we have used the same bound on 1^421 as in the previous section. Let 
Nq{£) C Nq denote those points {zi,Z2,zs) such that jzsl > e. The above 
calculation shows that gs{X) C -/Vo(^^). We compute the differential 



2k 

de 



1/Z3 -Zi/Z^ 



1/Z3 -Z2/zi . 

Using the fact that l^ip + |5;2p = l-^sPj we see that the sum of the squares of 
the matrix entries of d6 is 3/|z3p. It follows that 

Here V is any vector tangent to y. 

Therefore, given X G S^, the path gs{X) has speed at most ^ times 
the speed of the path gs{X), at corresponding points. An exercise in calculus 
shows that < .5 and 1^42 (s)| < .5 for all s G J. Prom this it follows 

that gs{X) has speed at most 1. Hence, the path gsiX) has speed at most 
This is a reformulation of Lemma 8.7. □ 

Let X G H U dQ^is, t)) be arbitrary. By Corollary 8.5, there is some point 
y G H U dQ,A{s) such that p{x,y) < a^. By Lemma 8.3, and the triangle 
inequality, the points x and y are both at least h — a> a from 5^ — E3. Hence, 
G F3. Prom Lemmas 8.6 and 8.7 

Prom the second statement of Lemma 8.2, and the triangle inequality: 
p{9t{x),gs{y)) < pigt{x),gs{x)) + p{gs{x),gs{y)) < 4aj + 2\s-t\ < b. 



9. The discreteness 

9.L Five disks. The set of disks 
rameter, just as in Section 7. Here we 
contains Wabc{s), for all s G J. This 
defined below. 



proof: General case 

WABcis) can be defined, for each pa- 
define a single collection WabCi which 
collection is the union of the five disks 
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Lemma/Definition. Let Wa{s) be the disk in W such that Ag c 
dWA{s). Then Wa{s) C Wa{s') if s > s'. In particular, Wa = Wa{s) is 
defined so that Wa{s) C Wa for all s E J. 

Proof. Here is the hyperbolic geodesic containing the fixed points of 
h-^-^s and which move together symmetricahy along dW2 as s increases. 

□ 

Lemma/Definition. Let W^{s) be the disk tangent to W2 at the fixed 
point p±^s of h±^s, Oj^d also tangent to dW . Let a = (s + s)/2. Let Wb be the 
smallest disk, concentric with W^{a) which contains both W^(s) and W^{s). 
Then /i±,s(W2) C for all s e J. 

Proof. By construction, /i±,s(W2) C W^{s). Since the disks W^{s) are 
all isometric they are all Euclidean isometric to each other, and move mono- 
tonically through the annulus Wi, we see that W^{s) C for all s e J. 



Lemma/Definition. Let W^{s) = h±^sh^.siW^). Then W^{s) c 
W^{s') ifs< s'. IfW^ = W^{s) then W${s) C W$ for all s e J. 

Proof. We consider the case of Wq, the other case following from sym- 
metry. Now, /i+,s/i_,s({0} X R^) is determined by its endpoints 



These arc the functions from Equation 23. An exercise in calculus show 
that the real parts (and hence the arguments) of the second coordinates of 
the above points are monotone for s £ J. Explicit evaluations at the end- 
points of J show that < arg/i+/i_(0, 1) < arg/i_|_/i_(0, — 1) < tt/2, that 
arg^+/t_(0, 1) decreases with the parameter and that arg/i+^_(0, — 1) in- 
creases. Hence Wq{s) C Wq{s) for all s E J. □ 

9.2. All words but two. Let a and b be as in the year lemma. 
Lemma 9.1 (Estimate 1). For all s G Jo, 

1. uwmA{s)r)nwABc = ^, 

2. (S U di},A{s)f C S3, the set from Section 8. 

Proof See Section 12. □ 



□ 
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Corollary 9.2. Ilw{^A{s,t)) n Wabc = 0- 

Proof. The same argument as in Corollary 8.5 implies that ilA(s,i) C 
(riA('S))"- The rest of the argument is immediate. □ 

The same symmetry as in Section 7 imples that Il^y t)) H Wabc = 0- 
As in Section 7, let Ui^t be the component of — Zi(s,t) which contains both 
fixed points of h-^-^th-,t- Let U2,t = d{Ui^t)- Once we establish the following 
lemma, the arguments in Section 7 apply verbatim to prove that h{U2,t) C Ui^t 
for all h ^ Hf — i^2,t- The reader should bear in mind, while reading our proof, 
that the thin set B{s,t) plays an important topological role in it. 

Lemma 9.3. Ui^t f^2,t- 

Proof. Assume this result is false. Let x be as Lemma 7.7. The same 
argument in Lemma 7.7, combined with Lemma 9.2, implies that there is a 
loop 7 G 5^ — r2(s, t) which links either dA or ddA, or both. Assume, without 
loss of generality, that it is the former. B{s,t) is a thin solid torus containing 
dA as its core curve. By construction, r] = dQ,'^{s,t) C B{s,t). This latter 
curve runs nearly parallel to OA, and is much longer than the cross sectional 
diameter of B{s,t). Hence, ij is nontrivial in the first homology of B(s,t). In 
particular, 7 and rj are linked. The same argument as in Lemma 7.7 implies 
that 7 intersects ^'^{s,t), which is a contradiction. □ 

9.3. Words of length 2. We found the following barrier function fg : 
^ H experimentally: 

(58) fs{u) = tis) + (.21 + .ll{s-s)) {z\s), 2t{s)) = QsiO- 

Here ( is as in Section 7. We choose the branch of ^ so that z{s) varies 
continuously. 

Let A(s) and Aj{s) be the sets associated to the barrier (fs,'^s), as in 
Section 7. Let Ii{s) = {— 00, t{s)] and his) = [t{s), 00). 

Lemma 9.4 (Estimate 2). Let c = a + b. For all s £ Jq and for all t G Ig, 

1. *,(S^U(di7A(s)") C Ai(s). 

2. 7rRO*,((dJ^d(s)r) C/i(s). 

3. ^s{Gs) C A2(s). Here Gs = {gs{E)r ^ {dsidQAisW. 

4. TTR o ^iJs{Hs) C his). Here H^ = {gs{d^d{s))y . 

5. TTR o '^stigtiPt)) > max fg. 



6. gt{pt)^{EUdnA{s))'^. 
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Proof. See Section 12. □ 

Corollary 9.5. Let s = Sk e Jo- For all t G Ik, 

1. ^sti^'JdnA{s,t)) c Ai(s). 

2. Trno'i!st{d^d{s,t)) C Ii{s). 

3. ^'stOfft(SUdJ^(s,t)) C A2(s). 

4. TTR o o gtid^ldis, t)) C /2(s). 

5. TTR o '^st{9t{Pt)) > max/5. 
6- gtiPt) 0SU(iriA(s,t). 

Proof. Statement 1 of the year lemma, combined with the triangle in- 
equality, implies that 

We also know that "^gt = o . Hence, the set defined on the left-hand side 
of Statement 1 of the corollary is contained in the set defined on the left-hand 
side of Statement 1 of Estimate 2. This proves Statement 1 of the corollary. 
Statements 2-4 of the corollary have very similar proofs. Statement 5 of the 
corollary has just been copied down from Statement 5 of Estimate 2, for the 
sake of exposition. Statement 6 of the corollary follows from Statement 6 of 
Estimate 2, and Statement 2 of the year lemma. □ 

Our goal is to show that Z(s,i) fl gt{Z{s,t)) = for i G 7^. We assume 
t < s. Once this is done, the same argument as in Section 7 shows that 
{Ui,t, U2,t) is a compressing pair for pt. 

Lemma 9.6. Z(s, t) n 5t(Z(s, t) c t) n gt{^'A{s, t)). 

Proof. By Lemma 5.3, 

- {pt}) C ^stid^^M) = ^s{d/S) C Ai(s). 
Combining this with Statement 1 of Corollary 9.5, we have 

*,t(Z(s,t) -{pt}) c Ai(s). 
Combining this with Statements 2 and 6 of Corollary 9.5 we have 

z(s,t)n5t(sudJ^A(s,t)) = 0. 
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The element r conjugates gt to g^^ and preserves every relevant set. By sym- 
metry, then, 

gt{Z{s,t))n{ElJdnA{s,t)) = 0. 
The lemma now follows from the definition of Z(s,t). □ 

Lemma 9.7. gti^^i^, t)) n n'^{s, t) = 0. 

Proof. First of all, we observe that gt{pt) / Pt, and Cl'^{s,t) D Et = {pt}- 
Also, gt{n'^{s,t)) nEt = {gt{pt)}- Hence, 

Pt^gt{n'd{s,t))nn'a{s,t). 

By Lemma 5.3 and Statement 2 of Corollary 9.5, 

TTno ^stWs,t) - {pt}) = TTno ^st{dn2s,t)) c TTRO ^I>,t{dnAis,t)) c Ii{s). 

Similarly, by Statement 4 of Corollary 9.5, 

TTjio^stO gt{dn'a{s,t)) C TTRO^^tO gt{dnA{s,t)) C his). 

Applying Statement 5 of Corollary 9.5 and Lemma 3.4 to each of the foliating 
arcs of gt(i}'^{s,t), we get 

7rB.o^,togt{n'a{s,t)) C his). 

Since /i(s) n his) = 0, the equations above combine to prove this lemma. □ 

Lemma 9.8. Let s E Jq, the set defined in Section 8. The hybrid disks 
(ihis) and dvis) are tame. 

Proof. We check the hypothesis of Lemma 5.7 for each parameter in Jq. 
(Compare Equation 43.) See Section 12 for details of this calculation. □ 

Lemma 9.9. ttr o ^'^ o gsiC) > max fs, for all s E Jo- 
Proof This is an explicit calculation. See Section 12 for details. □ 
Lemma 9.10. gti^'^is, t)) n fi^(s, t) = ^ forye {h, v}. 

Proof. Wc will take y = h. The case y = v has the same proof. We will 
show that the tame hybrid disks and the hybrid sectors 

satisfy the hypotheses of Lemma 5.8. 

1. The axis for ^'his,t) is iEt,pt). The axis for gti^'^is,t)) is iEt,gtipt))- 
Since t s, and since the fixed points of gt are not contained in the 
Clifford torus, we have qt = gtiPt) 7^ Pt- This shows hypothesis 1. 
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2. Since dQ,'^{s, t) C E(s, t) and dgt^d{s, t) C 5j(H(s, t)), Lemma 9.6 implies 
that 0,d{s,t) and gt{^h{s,t)) satisfy hypotheses 2 and 3. 

3. Note that ^^(C) C gs{d^ld{s))- Let a be the fohating arc of ys(drid(s)) 
which contains (. Lemma 9.9 imphes that ttr o ^ o gs{C) > max/^. 
Statement 5 of Estimate 2 imphes that ttr, 0*0^5(^5) > max/^. Lemma 

5.5 now imphes that ^s(a) C A2(s). From Lemma 5.3 and Estimate 2 
wc have ^,5(f7(_Bs,ps; Cs) — {ps}) C Ai(s). Since Ps a, we see that 
a n Q{Es,ps; Cs) = 0, which is hypothesis 4. |--| 

Lemma 9.1L gt{^x{s, t)) n Oy(s, t) = for x,y e {v, h}. 

Proof. The proof is the same as in the previous case. Ah we have to do 
is show that there is some fohating arc of gs{^y{s,t)) which is disjoint from 
Q,x{s,t). This foUows from Estimate 2 and Corohary 5.6. □ 



10. Numerical analysis 

10.1 The Clifford torus case. Recall that S is the part of the dented torus 
which is also part of the Clifford torus; S is the "outer rim" in Figure 10.1. 
We subdivide H into five Clifford rectangles, as shown in Figure 10.1. (What 
looks like twelve rectangles is really 5.) Let {Ej : j = 1,2,3,4,5} be the set 
of these pieces. Let E{t, u) be the rational parametrization of Ej, as described 
in Section 2.7. 

Ml Mn MT Ml 

^1 ^1 



£-4 



^4 



^5 



^5 



Hi HT Ht Mi 



Figure 10.1 

Given a rectangle Q C [0,1]^, let Ej{Q) be the Clifford subrectangle of 
Ej parametrized by Q. We now construct a ball B{Ej,Q) C such that 
Ej{Q)cB{Ej,Q). 

Let m be the midpoint of Q. Let V be the set of vertices of Q. Let B{Ej, Q) 
be the smallest ball centered at Ej{m) and containing the points Ej{V) in its 
boundary. All four points are contained in the boundary by symmetry. Since 
each C-arc foliating Ej is contained in a semicircle, we have Ej{Q) C B{Ej, Q). 
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The rectangle algorithm computes B['Rj,Q), returns a 1 if the t-curves of 
E]j((5) are longer than the it-curves and returns a 2 otherwise. The return of 
1 tells the computer that it can more efficiently decrease the size of (Q) by 
cutting Q in half so as to shorten the t-curves. A return of 2 tells the computer 
the opposite. 

Remark. The analysis above uses, in a crucial way, the fact that the cir- 
cular arcs in question are contained in semicircles. This property plays an 
important role in our analysis on several occasions. The key fact about such 
arcs is that their overall size is comparable to the distance between their end- 
points. This is not true for circular arcs in general. 

10.2. The hybrid sector case. Given an arc 7 C S^, let £'(7) denote the 
infimal e > such that every point of 7 is within e of either endpoint of 7. For 
instance, if 7 = [0, 1] C R then ^(7) = 1. 

Let CI = Q{E,p;S) be a hybrid sector, with parametrization fl(t,u). Let 
T = [ti, t2] and U = [ui, U2] be subintervals of [0, 1] . Define 

(59) n{T,u) = U n{t,u), n{t,u) = [j n{t,u), 

where ^{t, U) is an R-arc. We call ^{T, u) a transverse arc. The goal of this 
chapter is to produce computable bounds Lq {t, U) and (T, U) such that 

1. E{n{t,U)) < Ln{t,U). 

2. E{n{T,u)) < Wa{T, U) for any ueU. 

Our estimates are not general ones. We have only verified them for the clean 
subsectors of the hybrid sectors which actually participate in Estimates 1 and 
2. (We will define below what we mean by clean.) 

Let Ti and T2 be the two intervals obtained from splitting T in half. 
Likewise define Ui and 1/2. Let m = {ts,us) be the midpoint of Q. We have 

p{U{t, u),Cl{m) < p{0,{t, u),Cl{t3, u)) + p{Cl{t3, u), Q(t3, U3) 

< max(VF(ri, U), W{T2, U)) + max(W^(t3, Ui), W{t?,, U2)) = r{T, U). 

Thus, Q.{Q) is contained in the ball B{^l,Q) centered at Q,{m) and having 
radius r(T, U). 

The rectangle algorithm produces the ball B{Q, Q). Also, it returns a 1 if 

Ln{ti, U) < Wn{T, U), and returns a 2 if Ln(ti, U) > Wn{T, U). The integer 
values returned have a similar purpose as in the previous section. 

10.3. Distortion of stereographic projection. Define 

(60) U{e) = {{z,w)eS^\\z + l\<e}. 
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Note that 17(2) = S^. Let Bq be the map from Equation 7. For any smooth 
map T, let dT be the hnear differential of T. 

Lemma 10.1. If V is a unit tangent vector to 5^, based at a point p G 
U{e), then \\dBo{V)\\ > l/{eV3). IfV is £-integral then \\d{Trc o Bo){V)\ > 
1/e. 



iw{z + 1) 



Proof. By rotational symmetry, continuity, and integration, it suffiees to 
assume that Bo(p) € R x R. Note that Bg ^(R x R) consists of those points 
{z, w) such that w/ (z + 1) G R. 

For {z,w) e Bo^(R* x R), define 

{1+zf \-W] \ " 

"""-^^TVi z J' '''=[i{z + l)_ 

A vector u is tangent (resp. £^-integral) to at p if and only if 3f?(p, u) = 
(resp (p,!^) = 0.) Using this criterion, we see that all three vectors above are 
tangent to S^, and that is also iS-integral. 

Observe that all vectors have norm at most e. In the third case, this is 
a short exercise in calculus, which uses the fact that 12:^ + 1^1^ = 1. We will 
show that dBo maps the vectors above to an orthonormal basis. This implies 
our first differential bound. We will see, in particular, that d(7rc o Bq) maps 
Vr to a unit vector. This implies the second bound. 

Let Pi and (32 be the component functions of Bq. We compute 

1 



d(3i = 
dPiiue) = 

dhii^t) = 



—w 

(1 + Z)2'T + ^ 

i; (i/3i(t't) = 0, 



± lim — 9= 

e^O e 



1. 



dPi{vr) = 1; 

"1-z l-z-ie(z + 1)2/2" 
l + z ~ l + z + ie(z + 1)2/2. 

As we said above, vg is a real multiple of (0, iw), and this latter vector generates 
the flow (-2, ty) — > (z,exp(ii) w). The map Bq conjugates this flow to the flow 
on Ti generated by r^. Thus, we see that df32{^e) = 0- The ^^-integrality of 
d(3i{ur), together with the fact that d(3i{ur) = 1, implies that df32{i'r) = 0- D 

Lemma 10.2. Suppose {z,w) € S'^ — U{e). If V is a unit tangent vector 
based at {z,w) then \\dBo{V)\\ < Ve^^ ^ 5^-4 < ^ 

max(l,e ). 

Proof. Let Bi be such that Bq = vr o Bi. Here Tr{z,w) = {z,Q{w)). Note 
that ||dBo(y)|| < ||ciBi(F)||. We compute 

'l/il + z) w/{l + zf^ 
.2/(1 + 2)2 



dBi 



The lemma follows straightaway. 



□ 
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10.4. Clean hybrid sectors. Let /n be the endpoint map for ^l{E,p] S), as 
defined in Section 4.1. We say that Cl{E,p; S) is clean if 

1. S is contained in a semicircle, 

2. fn{S) C -B is contained in a semicircle, 

3. fci is injective on S. 

Let y be a Heisenberg stereographic projection mapping Q, to standard 
position. Jq is injective on S if and only if VTRoy is injective on S. Now, Y{S) 
is contained in an ellipse, so that ttr has at most one local maximum and one 
local minumum on Y{S). 

The cleaning algorithm does the following: 

1. Finds the points {pi} of S where /q is extremized. 

2. Sorts the points {m} U {pi}, according to the order they appear on S. 
Here m is the midpoint of S. 

3. Subdivides S into subarcs {Si}, bounded by the sorted points. 

4. Checks that the hybrid sectors {i},{E,p; Si} are clean. 

We implement the first step using a simple iterative method. We will 
explain how the method finds the point /v, where ttr o Y attains its minimum. 
Finding the maximum point is similar. We first locate n roughly, and thereby 
produce a map (pQ : [0, 1] — ^ 5 whose image contains /x. Next, we evaluate the 
points ttr, o y o (a(j/3)) for j = 0, 1, 2, 3. The minimum occurs either at 1/3 
or 2/3. In the first case, we define = (j)o o Ai. Here Ai is the affine map 
which takes [0, 1] to [0,2/3]. In the other case, we define = i^o o ^2- Here 
A2 is the affine map which takes [0, 1] to [1/3, 1]. Iterating, we get a sequence 
of maps (/)o, 01, (/>2"- whose images rapidly converge to //. 

Let {Si} be the set of arcs produced by the first three steps of the cleaning 
algorithm. We will implement step 4 of the cleaning algorithm by showing p 
is farther from either of the two points of fnidSi) than these points are from 
each other. 

10.5. Estimating the R-arcs. Let 0,{E,p;S) be a hybrid sector. Let St 
be the canonical parametrization for S. Let At = fn{St) be the image of the 
endpoint map, as in Section 5.3. We define the endpoint curves: 



We have used the group law on S^ to define these curves. More generally, we 
define 



(61) 



Ep{t) = Bo{-p-At'y, 



Es{t) = Bo{-SfAt'). 



(62) 



En{t,u) = -n{t,u)-A, 



-1 
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Lemma 10.3. En{t, u) = uEp{t) + (1 - u)Es{t). 

Proof. Apply Lemma 2.1 to the three points X = —p ■ A^^ and Y = 
-S ■ A^^ and Z = -At ■ A^^ = (-1,0). □ 

Combining Lemmas 10.1 and 10.3 we see that, for all u E U, we have 
p{^{t, u),^l{t, ui) < Li{t, U), where 

(63) Li{t,U) = [2] \ui-U2\ \7rc{Ep{t)) - irciEsm- 

The 2 above is bracketed for purposes which will become clear in Section 10.7. 

10.6. Estimating the transverse arcs. A dyadic interval is an interval of 
the form 6/2^], for a, b, k integers. When we run our main algorithm in 

the next section, we will only have to deal with dyadic subintervals of [0, 1]. 
Our analysis of transverse arcs works for dyadic subintervals which have length 
at most 1/4. If T is a longer interval, we set W{T, U) = tt, the diameter of S^. 
Henceforth, we assume that T is a dyadic subinterval of [0, 1], having length 
at most 1/4. 




Figure 10.6 

We would like to compare = ri(i3, u) and uj^ = Q,{t4,u), for ta, t4 G T. 
The point ooj lies in the R-arc Qj = Q,{tj, [0, 1]). Let Aj = At- and Sj = Sty 

Let Rj be right multiplication by —A~^. Let = Rj{Xj) for any point 
Xj G Qj. In particular, let = Rj{0,j). We have A[ = A2 = (-1,0). Since 
Rj moves all points on S"^ the same distance, and since Ai, A2, A3, A4 are 
contained, in order, on a semicircle, we get 

(64) p{uJ3,UJ4) < P(^^3,^4) + P{A3,A4) < p{uj'3,Uj'i) +^(^1,^2)- 

Let Bq be the stereographic projection from Equation 7. Note that 
Bo(0,-l) = GO. Let X" = Bo(X') for all relevant X. We equip H with 
the Euclidean norm. In subsection 10.8 we will prove the following result. 
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Lemma 10.4. If is a dean hybrid sector, produced when the cleaning 
algorithm is applied to a hybrid sector participating in Estimates 1 and 2, then 

\\p'i-p'i\\ < .002 + \\p';-p% \\S'^ - S'lW < .002 + \\S'I - 



By Lemma 2.1, uJj = up'- + (1 — u)S'-. Combining this equation with 
Lemma 10.1, and Lemma 10.4, we see that p{u''^,oj'^ < Wi{T,u), where 

(65) Wi{T,u) = [2]V3(.002 + tx||Bo(i?i(pi))-Bo(i?2(P2))|| 

+ (1 - u)\\BoiRi{Si)) - BoiR2{S2), 



Again, 2 above is dehberately bracketed. The coefficients in the last equation 
are monotone for u G [0, 1]. Hence, independent of the choice oi u E U, we get 
p(a;3,wl) < Wi{T, U), where 

(66) Wi{T,U) = ma^{Wi{T,ui),Wi{T,U2)). 



10.7. Bootstrap argument. Using the bounds above, we produce a ball 
B'{Q,,Q) such that n{Q) C B'{Cl,Q). Our calculation of the radius of B' 
depends only on the fact that fl{Q) C ^7(2). Given B' , wc find the best £ > 
such that B' C U (e). Then, using the full force of Lemma 10.1, we can replace 
each occurrence of the bracketed 2 above by e. The smaller ball B{yi,Q) 
also contains ^{Q). This bootstrap argument reduces our computation time 
roughly by a factor of 10. 

10.8. Proof of Lemma 10.4. Let 7 : [0, 1] ^ H be a curve. Given a dyadic 
subinterval / C [0, 1] let .6(7, /) be the ball in 7i which has, as a diameter, the 
segment joining the two endpoints 7(<9/). Let Im be the midpoint of I. Say 
that / is good if j{Im) C -6(7, /). 

Let T denote the set of dyadic subintcrvals of T. Say that a Lipschitz 
estimator for 7 is map A : R such that ||7(i) — 7(/m)|| < A(/) for all t & I. 
Consider this algorithm: 

1. Let LIST be the partition of T consisting of just T. 

2. If LIST is empty, halt. Otherwise... 

3. Let / be the first interval on LIST. If A(/) < e, delete / from LIST and 
go to Step 2. Otherwise... If / is good, replace / by the two intervals 
obtained by cutting / in half, and go to step 2. Otherwise fail. 
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If this algorithm halts without failure, the information constitutes a proof 
that 7(t) is contained within the e neighborhood of B{'j,T). The point is that 
there is a sequence T = To,Ti, . . . ,Tk of dyadic intervals such that \\j{t) — 
jidTk)\\ < e and B{^,Tk) C 5(7,Tfc_i)... C B{^,Tq). 

When we run this algorithm, with the Lipschitz estimators defined below, 
fore = .002, and for the curves (i) = ^Q{-p-A^^) andS'"(t) = Bo(-S'rA~^)' 
it halts without failure for every parameter in Jq\J{s}. See Section 12 for details 
of the calculation. This establishes the lemma. 

Remark. For several parameters s G Jq, one of the hybrid sectors pro- 
duced by the cleaning algorithm is tiny. In these cases, checking the "control 
by balls" condition involves the comparison of two tiny quantities. The round- 
off error introduced by the computer interferes with this comparison. In these 
exceptional cases, we simply use our Lipschitz estimators to verify, in advance, 
that the relevant curves are entirely contained in a ball of radius .002, thereby 
bypassing our algorithm. 

Lemma 10.5. Let m be the midpoint of [a,b]. Let 
Ai = mayi{p{Aa, Am), p{Ab, Am)); M = |1 + 7ri(-p • ^~^)| -Ai. 

The quantity Ap([a,5]) = Ai\/6max(l, A2 ^) is a Lipschitz estimator for p" . 

Proof. Let m be the midpoint of [a, 6]. Since S\t is contained in a semi- 
circle, ) < Ai. By the triangle inequality p{p'^,Pm) < Ai. By 
definition, = 'Bq^{p'^) and likewise for c. Therefore, 

p{B^\p"{c)),B^\p"{m))<X,. 

By the triangle inequality, 

|7ri(Bo Hp" (*))) + 1| > A2; te[c,m]. 

The points BoHp"(a)) and B^^{p"{c)) can be connected by an arc having 
length at most Ai. By Lemma 10.2, the restriction of Bq to this arc expands 
distances along this arc by at most -\/6max(l, A^^). □ 

Essentially the same argument gives 

Lemma 10.6. Let m be the midpoint of [a,b]. Let 

Ai = mayi{p{Aa,Am),p{Ab,Am)) + iaayi{p{Sa,Sm),p{Sb,Sm)), 

A2 = \i + M-s-A-^)\-Xi. 

Then Ap([a,6]) = Aiv^max(l, A2 ^) is a Lipschitz estimator for S" . 
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11. Verifying the estimates 

All our estimates work, one parameter at a time. Suppose the parameter 
is fixed. The first step, in both cases, is to split the relevant hybrid sectors into 
clean hybrid sectors, and (for Estimate 2, Lemma 9.4) to split H into {Hj}, as 
in Section 10.1. 

Henceforth, we will say that a piece is either a clean hybrid sector par- 
ticipating in Estimates 1 (Lemma 9.1) and 2, or one of the Hj. This chapter 
discusses how we verify that Estimates 1 and 2 hold for an individual piece. 
We will also explain the remaining statements in Estimate 2, which do not 
pertain to individual pieces. 

11.1. The subdivision algorithm. Let Qq = [0,1]^. Given a rectangle 
Q = [ai, hi] X [a2, 62] C Qo, let 

(67) Si(Q) = [ai,(ai + 6i)/2] X [61,62] U[(ai + 6i)/2, 02] X [61,62], 

(68) S2(Q) = [ai, 02] X [61, (61 + 62)72] U [ai, 02] x [(61 + 62)/, 62]. 

These are the two dyadic subdivisions of Q into two rectangles. If L = 
{Qk} is a finite list of rectangles, let T,j{L) be the list obtained from L by 
deleting the last rectangle Qn and appending Ilj{Qn). 

Let Q denote the set of rectangles of [0, 1]^. Say that a subdivision test is a 
map T : {0, 1, 2}. We are interested in proving the existence of a partition 
of Qo into rectangles {Qi, ■ ■ ■ , Qn} such that T{Qj) = for j = 1, . . . , n. Here 
is the subdivision algorithm: 

1. Let LIST be a list with Qo as its only member. 

2. If LIST is empty, halt. If not, let Q be the last rectangle in LIST. 

3. If T{Q) = delete Q from LIST and go to step 2. 

4. If T(Q) =je {1, 2} then replace LIST by Sj(LIST) and go to step 2. 

If the subdivision algorithm halts, the information constitutes a proof that the 
desired partition exists. 

11.2. Estimate 1. Let a and 6 be as in Estimate 1. Let be a piece. Let 
Q C [0, 1]^ be a rectangle. Let Pq and Sj. be as in Equation 56. Here is the 
subdivision test: 

1. Let B = B{p,,Q) be the ball produced by the rectangle algorithm. Let 
(c, r) be the center and radius of B. Let j G {1, 2} be the integer produced 
by the rectangle algorithm. 
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2. For i = 1, 2, 3, 4, 5, let {(q , r^)} be the collections of centers and radii for 
the disks in Wabc- If mini(|nt^(c) — Ci\ — r — ri) < a, let T{Q) = j. 
Otherwise... 

3. If (||nvi/od(c))-Po(s)||-|-r < 6, letr(Q) = j. Otherwise let r(Q) = 0. 

Po{s) is evaluated from Equation 56 (§8.3). See Section 12.6 for a list of 
the centers and radii of the disks used in this estimate. 



11.3. Variation of the barrier function. We use the notation of Section 9.3. 
Let fs be the barrier function in Equation 58. Let tti : jS^ x R ^ S*^ be 
projection. Let ^'^ be the elevation map associated to the parameter s. For 
q £ S'^ define 

(69) Fs{q) = /s0 7rio*,(g)+7rRO^',(g), 

(70) Gs{q) = 7rRO*,(g)-i(s). 

We have defined this function so that 

1. ^sip) G Ai(s) if and only if Fs{p) < 0. 

2. ^sip) G Mis) if and only if Fs{p) > 0. 

3. TTR o ^sip) G h if and only if Gs{p) < 0. 

4. TTR o *s(p) G h if and only if Gs{p) > 0. 

Let B C he a ball with center c and radius r. We define Fs{B) = 
[j^^s Fs{b). If p G S we redefine Fs{B) = [—10,10], which causes B to fail 

all further tests applied to it. If c = 0, we replace B by the ball of radius 
r + .0000001 centered at .0000001. Henceforth, we assume that c 7^ and 
p ^ B. We also assume that the parameter s is fixed. Let c be the affinely 
normalized lift of c. 

Given a vector V = (zi, 22, 23), let 



(71) \\Vh = ^J\zi\^ + \z2\^. 
Define 

(72) L+{B,V) = \{c,V)\+r\\V\\2; L_{B,V) = max{0,\{c,V)\-r\\V\\2). 

Let ps and qg be any lifts of the points Ps,qs € Eg. Let E* be any vector 
polar to Es- Let r be the radius of B. Define 

(73) a, = minf2vr,^^M^+ ^I^^^H^ 



L^{B,E*)\ \L-{B,ps 

(74) £1 = min(ai,2)/3(s), /3(s) = .21 + .ll(s - s) 
^ A\Qs\\2 , r||p^||2-L+(-B,g^) 

(75) £2 = r /p „ \ + 



L^{B,ps) L^{B,psf 
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Lemma 11.1 (the barrier test). With the above notation, 

1. If Fs{c) + ei + £2 < then ^ s{B) d Ki. 

2. IfFslc) - £1 - £2 > then ^{B) c A2. 

3. // Gs{c) + £2 < then ttr o ^s{B) C h. 

4. // Gs{c) - £2 > then ttr o ^ s{B) C h- 

5. //Gs(c) +£2 + /3(s) < then s{B) C Ki. 

6. //G^jc) - £2 - /?(s) > then s{B) C A2. 



Proof. Let g G i? be an arbitrary point. We define {zq,tq) = *s(g). 
Likewise {zc,tc) = ^s(c). We first prove that 



(76) 



A = I arg(zc) - arg(2;g)| < ai. 



To see this, let 7 : [0, d] he a constant-speed geodesic which connects c 

to p. Let 7 be the affinely normahzed Hft of 7. We certainly have A < 27r. 
Define X{t) = arg(5' o 7(i)). By integration, 



A < sup ||X'(t)||5 < sup \\X'{t)\\r. 
t t 



By definition, 

X{t) = axg 



QL{t); L{t) = log 



Here log is a branch of the complex logarithm. By calculus. 



\X'{t)\ < \\L'{t)\\ 



+ 



m),Et) 

The third coordinate of 7'(t) is zero and ||7'(i)|| = 1- By the Cauchy-Schwarz 
inequality, 

|(7'(i),^)l<mi2. 

Here V is any vector whose third coordinate is 0. Another application of the 
Cauchy-Schwarz inequality gives 



Therefore 



\{^{t),V)\>L_{B,V). 



II ^* II 

\x'(t)\< _ ii.J'i . + 



L^{B,E*) L^{B,p) 



£1. 



Equation 76 follows straight from this. 
A very similar argument proves 



(77) 



tc\ < £2- 
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Given the definition of we have 

(78) max/, -min/, = 2/3(s); \dfs/du\ < P{s). 
The first inequahty in Equation 78 immediately gives 

(79) /s o vri o ^,(c) - /, o TTi o q>,{q) < 2/3(s). 

Integrating the second inequality in Equation 78 over the portion of which 
contains tti o "i/siB), we get 

(80) /s o vri o ^'.(c) - /, o TTi o ^,{q) < ai/3(s). 

Combining Equations 77, 79 and 80 we see that \Fs{q) — -Fs(c)| < ei + £2- 
Hence Fs{q) < 0. This is equivalent to the statement that ^'(g) C Ai. The 
other statements have similar proofs. □ 

11.4. Acting on the ball. Let and Po{s), and i^^ be as in Lemma 8.6. 
Let B & he a ball, with radius r and center c. 
Define 

1 , 2k 



(81) k = , K = min(4,— V27 + 31fc2). 

Fi(c) — -P(s)| — r 27 



Lemma 11.2. gs{B) c X{B,gs), where X{B,gs) is the ball having center 
gs{c) and radius Kr. 

Proof. By construction, every point x & B has the property that 
(c, x) G -Ffe. Therefore, the Lipschitz constant of g\B is at most K by Lemma 8.6. 

□ 

11.5. Confining a particular curve. Let s G Jq and let be the corre- 
sponding partition interval. Let Ps be the point used in the definition of the 
hybrid sectors of Vt/^{s). 

Lemma 11.3. Let (3s be the ball of radius .001 centered at gs{Ps)- Then 

{gt{Pt)T e (3s for all t e is- 

Proof. Recall that a = .00046. Let Po(s) be as in Equation 56. Define 

(82) 

k= , . , ,1 , K, = -V27 + 31k^; K2 = -\Is\. 

|7ri(ps) - Po{s)\ - 2a - .02 27 3 

From the year lemma, p{ps,pt) < a. An exercise in calculus shows that the 
real and imaginary parts of Pois) are monotone in s. Evaluating at the end- 
points, we see that \Po{s) — Po{t)\ < .02 for all s,t G J. By the triangle 
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inequality, \7ri{pt) — -Po(*)| < {T^iiPs) — Po{s)\ + -02 + a. We have chosen our 
constant k so that the geodesic in which joins Ps to pt remains in Sj^.. Hence, 
iPs,Pt) G Ft. Applying Lemma 8.6 we see that p{gt{Ps), 9t{pt)) < Kia. Apply- 
ing Lemma 8.7 we see that p{gsiPs), gt{Ps)) < By the triangle inequality, 
Pidsips), 9t{Pt)) < -f^ifl + -f^2- An explicit calculation, for all s G J shows that 
the right-hand side of this equation is at most .00054 = .001 — a. See Section 12 
for details. □ 

Corollary 11.4. Statement 6 of Estimate 2 {from Lemma 9.4) is true. 

Proof. Let m+(s) be the maximum value on attained by fg. Let 

Xs = Yi\0,m+is)). 

Here Yg is the Heisenberg stereographic projection used to define ^g- Let 
Tt = ttr, o "^gt ° 9t- We compute explicitly that Tg{pg) > max/ for all s G 
Jq. (See §12 for details of the calculation.) We know that gt{pt) 7^ Pt, and 
Tt{pt) = oo- Thus, Tt{pt) 7^ oo. To verify Statement 6, we just have to rule 
out the possitility that Tt{pt) = m+(s). If Tt{pt) = m+(s) then Ygt o gt{pt) = 
YgO<j)~togt{pt) = Yg{xg) = (0,m+(s)). In other words, (t)gt{xs) = gt{pt)- Hence, 
Xg E Pg. For each parameter s G Jq we compute explicitly that Xg ^ j3g. See 
Section 12 for details. □ 

11.6. Estimate 2. Let $7 be a piece. Let Q C [0, 1] be a rectangle. Here is 
the subdivision test for Estimate 2. 

1. Let B = dB{Q,, Q), where B{Q,, Q) is the ball returned from the rectangle 
algorithm. Let j G {1,2} be the integer returned from the rectangle 
algorithm. 

2. (Statement 1) Use the Barrier Tests to see if ^{B"') C Ai. If not, set 
T{Q) = j. Otherwise... 

3. (Statement 2, when relevant) Use the Barrier Tests to see if ttr, o \I'(iJ") 
C h. If not, set T{Q) = j. Otherwise... 

4. Let Bg = X{gg,B). 

5. (Statement 3) Use the Barrier Tests to see if ^{Bg) C A2. If not, let 
T{Q)=j. Otherwise... 

6. (Statement 4, when relevant) Use the Barrier Tests to see if ttr, o ^(Bg) 
C I2. If not, let T{Q) = j. Otherwise... 

7. (Statement 5) If B n /?s 7^ set T{Q) = j. Otherwise... 



8. Set T{Q) = 0. 
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12. Implementation and roundoff error 



12.1. Overview. The bulk of our computing is done in C. We use C 
code for Lemma 10.4 and also for Estimates 1 and 2. The C program only 
uses mathematical operations which are governed by the IEEE standards [I]. 
Namely, 



2. The < and > operations; 

3. decimal <^ binary conversion. 

We implement our code using interval arithmetic, so that the computations 
themselves produce the error bounds. 

There are a number of places in the paper where we need a small amount 
of fairly precise computation. We perform these computations in Mathematica 
[W], which has arbitrary precision arithmetic. We compute all quantities to 
roughly 20 decimal places of precision. While Mathematica is not guaranteed to 
be free from computational bugs (and bugs have been found!) we think there is 
general agreement that Mathematica does not have bugs in the basic operations 
we use— the ordinary arithmetic operations and the numerical extraction of 
cube roots. We have ruled out such computational bugs by random tests on 
our data— e.g. plugging the root of a cubic back into the cubic. 

The C program relics on some initial values, depending on the parameter s. 
For instance, the C program needs to have the vector E* which in turn depends 
on the eigenvalues of gg. Our Mathematica code computes all initial values for 
the C code, and stores them in auxiliary files. When the C code runs, it reads 
in these files as needed. Thus, the C code never has to compute any quantity 
which depends on functions which are not listed above. 

We will first discuss certain features of our Mathematica code. Following 
this, we will explain how the computer represents real numbers, and how we 
implement interval arithmetic. At the end of this chapter, we give a record of 
all the calculations. 

12.2. The Mathematica code. Let Qg be the matrix in Equation 23. Let 
E* = (e(s), e(s), 1) be the affinely normalized eigenvector associated to Ag. 
Equation 30 gives e(s) as a linear function of A^. Computing e(s) accurately 
boils down to computing A^ accurately. Once we have an accurate value of A^, 
all our Mathematica calculations are straightforward implementations of what 
is explained in the paper. 

Finding the roots of the charactistic polynomial for gg leads to the formula 



1 



The +, — , X , and ^ operations; 



(83) 



A 



a{s) + b{s) + c{s) 
12(1 + s2) 
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where 

(84) a(s) = -68 + 64?s-4s^ 

8exp(47rV3)(-775 + 173Ais + 1203s^ - 308zs^ - ^Qs"^ 
+ 6is^ + / + 6^/3(l + 5^)3/ V3s3 - 125)^/3^ 
5440 - 7936ZS - 2688^2 + 256is^ + 64s^ 

This formula is vaHd, in particular, for s G J. Once Mathematica chooses 
branches of the square and cube root functions, it can compute the expression 
above to, say, 1000 decimal places. We will stop at 20 decimal places. 

Mathematica can consistently choose a branch of the square root function 
for all s G J. Mathematica takes the positive branch of Vl + and the 
positive imaginary branch of \/3s^ — 125. Mathematica needs to take a cube 
root in order to evaluate b{s). We have Mathematica compute b{s) along with 
Xg. For successive computed parameters si,S2 G J, Mathematica checks 

(87) |6(s,)|>50; \b{si) - b{s2)\ < 1. 

These bounds show that Mathematica cannot switch branches of the cube 
root when computing successive parameters. An explicit check shows that 
Mathematica chooses the correct branches for s. Hence, the correct branches 
are always taken during the computation. 

Once again, we ruled out computational bugs in the Mathematica code by 
randomly checking that the computed eigenvalues and eigenvectors are indeed 
eigenvalues and eigenvectors. 

12.3. Performing the cleaning algorithm. For each of the 366 parameters 
s £ Jq U {s}, we perform the cleaning algorithm on the hybrid sectors in 
^Ais). Each of the resulting clean hybrid sectors fl{Es,Ps; S^) is specified by 
the parameter s, the axis, which is the collection {E*,ps,qs,C} normalized as 
in Section 5.4, and an equalized choice of lifts {S^, S2, S^}. Here and S2 
are the endpoints of and is the midpoint of — . 

The axis is the same for all the hybrid sectors. The real and imaginary 
parts of the components of all relevant vectors are what is stored in the file. 
Everything is stored to 20 significant digits. (This is slightly more precision 
than is used in our C code.) The routine Clean [data.s , s] performs the clean- 
ing algorithm for the parameter s and stores the result in the file datas. The 
data contained in the created files is read in by the C code at the time of its 
operation. 



(85) b{s) = 

(86) c{s) = 
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12.4. Doubles. Our C code represents real numbers by (iou6/es. According 
to [I, §3.2.2], a double [x] is an object of the form (s,e, /). Here 

1. s is a single bit, determining the sign of x. 

2. e is an 11 bit word, representing an integer between and 2047. 

3. f = -bi - ■ ■ 652 is a 52 bit word. 
The real number represented by [x] is 



(The exceptions to this formula, which are detailed in [I, §3.2.2], do not arise 
in our calculations.) For example, 3f f 6a09e667f 3bcd represents the double 
closest to -\/23 The first bit of 3=0011 is 0. Thus s = 0. The concatenation of 
Oil and ff=llllllll is 01111111111. Thus e = 1023, and (_i)s2^-i°23 ^ i 
The word 6a. . . expands as 01101010 ... . Thus 



To access the words comprising a double, we introduce an exact: 

typedef struct 
{int a,b;} exact. 

Given an exact E, the integer E . a represents the first 32 bit word of the double 
(the high bits) and the integer E.b represents the second 32 bit word of the 
double (the low bits). Here is the code which enables one to convert between 
doubles and exacts. 

exact double-to-exact (x) 
double x; 



exact E; 

double *pointerl ; 
int *pointer2; 
pointer l=&;x; 

pointer2=( (int*) pointerl) ; 
E.a=*pointer2; 
++pointer2 ; 
E.b=*pointer2; 
return (E) ; 



(88) 





{ 



} 
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double exact-to-double (E) 
exact E; 

{ 

double x; 
double *poiiiterl; 
int *pointer2; 
pointer l=&;x; 

pointer2=( (int*) pointerl) ; 
*pointer2=E. a; 
++pointer2 ; 
*pointer2=E.b; 
return (x) ; 

} 



12.5. Interval arithmetic. Let D be the set of doubles. Define 

Ro = {x G R| |x| < 2^°23|^ 

Summarizing [I, §3.2.2, 4.1, 5.6], we see that there is a map Rq D, which 
maps each point x € Rq to some [x] G D which is closest to x. In case there 
are several equally close choices, the computer chooses one, as detailed in 
[I, §4.1]. 

Regarding the five basic operations, [I, §5] states that each of the opera- 
tions shall be performed as if it first produced an intermediate result correct to 
infinite precision and with unbounded range, and then coerced this intermediate 
result to fit into the destination's format. Thus, 

(89) ^/\x\ ■■= I^/HM)] [x]*[y]:=[ri[x])*ri[y])], *g{+,-,x,-}. 

We order exacts lexicographically, simply by concatenating the two inte- 
gers, treating them as a single integer. The inclusion r : D — R induces a 
linear order on the doubles. Beautifully, the linear order on positive doubles 
coincides with the linear order we have imposed on the corresponding exacts. 
Every double x (except the largest and smallest) has a unique successor x+ > x 
and a unique predecessor x_ < x. To compute these, we convert x to an exact, 
concatenate the integers, add 1 or —1, split apart the resulting integer, and 
convert back. 

One exceptional case requires mention: To avoid a certain kind of under- 
flow error, we redefined 0± to be the double closest to ztlO"^'^. In hindsight, 
this redefinition seems unnecessary. At any rate, it is certainly harmless to our 
calculations. 
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An interval is a pair / = (x, y) of doubles, such that x < y. Say that I 
bounds z G Rq if x < [z] < y. This is true if and only ii x < z < y. Define 



Let dil and 82! be the endpoints of an interval /. For j = 1,2, let Ij = 
{[xj], [yj]) be intervals. For * G {+,— , x,^} we define 



In all cases, one performs the relevant computations on the endpoints of the 
interval and then pushes the interval one click outward, to guarantee that 

1. If Ij bounds Xj then Ii * I2 bounds xi* X2- 

2. If / > bounds x > then s/l bounds s/x. 

The only exception to this rule occurs when we are computing I1/I2, and the 
two endpoints of I2 have different signs. In this case we automatically fail 
whatever computational test we are working on. 

A potential problem, related to the exceptional case, is that we might 
divide by a very small number, producing an overflow error. All the tests 
performed by our code fail if a calculation produces a sufficiently large number. 
Thus, an overflow error does not occur in the calculations which are relevant 
to our proof. 

We define a complex interval to be an expression of the form X+iY, where 
X and Y are intervals. We define a vector interval to be a triple of complex 

intervals. We define a matrix interval to be a triple of vector intervals. And 
so on. We say that a vector interval bounds a vector if the components of the 
vector interval bound the components of the vector. Likewise for the other 
structures. 

The algebra of the interval structures is identical to the corresponding 
algebra for the usual structures. At every step of our computation, the actual 
object is bounded by the corresponding interval version of the object. Thus, 
if one of our algorithms halts with success, as implemented, the information 
constitutes a proof that a perfectly accurate computing machine would also 
halt with success. 

12.6. Record of the calculations. Our code is contained in the directory 
res/Computers/Proof. This directory is contained in the University of Mary- 
land College Park mathnet system. 

The calculations from Section 3.4 are called CALCl and CALC2, and 
appear in the Mathematica code in the same order they appear in Section 3.4. 



(90) 



[x,y]o = [x-,y+]. 



(91) Ii* I2 = (rninSjJi * djL^ 
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The calculation CALC3 is the calculation for the year lemma. The calculation 
CALC4 establishes Lemma 9.8. The calculation CALC5 establishes Lemma 
9.9 and CALC6 establishes Lemma 11.3. Also, CALC7 establishes Corollary 
11.4. We performed the first three calculations on Sept. 2, 1999. The remaining 
four were done on Sept. 30, 1999. For all calculations, we used the computer 
legendre.umd.edu, which is a Sparc Ultra 5. 

Our Mathematica code also calculates the five disks used in estimate 1. 
We computed these disks on Sept. 2, 1999, on the same machine. Here is the 
result: 

1. Wa has center 1.0757057484009542760... and 
radius .3964124835860459412.... 

2. has center .8425731633320585661... ± (.13646924571449026840.. .)z 
and radius .1513081547068024.... 

3. has center 0.9316949906249123735... ± (0.3892740451354294927... )z 
and radius .13996370162172745806... 

We performed the C calculations, during the period from Sept. 24, 1999 to 
Sept. 28, 1999. We used, in parallel, the following Sparc Ultra 5 workstations: 

legendre . umd . edu 
noether . umd . edu 
kummer . umd . edu 
galois .umd. edu 
pascal.umd.edu 
sylow.umd. edu 
descartes . umd . edu 
maclaurin . umd . edu 
monge . umd . edu 
poisson . umd . edu 
Simpson . umd . edu 

These ghosts of famous departed mathematicians all took part in establishing 
our theorem. 

13. Computer plots 

Figure 13.1 shows Ilw{d^}ci{s))- The five grey disks are Wabc- The curves 
in the picture are individual foliating R-arcs, which attach to dW2- 

Figure 13.2 shows a close up of Figure 13.1, with some additional points 
plotted. 
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Figure 13.2 
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Figure 13.3 shows Iliy The five grey disks are Wabc- The curves 

in the picture are individual fohating R-arcs, which attach to dW2- 




Figure 13.3 



Figure 13.4 shows a close up of Figure 13.3, with some additional points 
plotted. 




Figure 13.4 
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Figure 13.5 shows 'n.w{dQ,fi{s)). The five grey disks are Wabc- The curves 
in the picture are individual fohating R-arcs, which attach to dW2- 




Figure 13.5 

Figure 13.6 shows a close up of Figure 13.6, with some additional points 
plotted. 




Figure 13.6 

Figures 13.7 shows ^-s{dQ,d{s)) and gs{^-s{d^(i{'s)))- The cylinder T has 
been identified to R/27rZ x Z via {u,t) (argM,t). The sine curve is the 
common boundary of Ai(/) and A2(/), with / as in Equation 50. 
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Figure 13.8 shows "i>s{d^ly{s)) and gs{"^s{d^v{s)))- When h replaces v, 
the picture looks exactly like Figure 13.8, but rotated 180° about the points 
of symmetry of A(s). 




Figure 13.8 
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In Section 10 we break H into 5 pieces: H = Hi U...UH5. (Compare Figure 
10.1.) Figure 13.10 shows ^'^(Ss U S3) and ^'5(5^(^2 U S3)). 




Figure 13.10 



Figure 13.11 shows a closeup of Figure 13.10. A somewhat different set of 
points has been plotted. 




Figure 13.11 
Figure 13.12 shows ^'s(Si) and ^'^(^^(Si)). 
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Figure 13.12 

We have r(H2 U S3) = H4 U H5. Thus, the picture for ^ji'^i U ^^5) and 
*s(5^(("4 U *5)) is obtained by rotating Figures 13.10 and 13.11 about the 
point of symmetry on A-j 

University of Maryland, College Park, MD 
E-mail address: res@math.umd.edu 
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